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Abstract 

In ongoing and upcoming hadron collider experiments, top quark physics will play an important role 
in testing the Standard Model and its possible extensions. In this work we present analytic results for 
the differential cross sections of top quark pair production in hadronic collisions at next-to-leading order 
^ in the QCD coupling, keeping the full dependence on the spins of the top quarks. These results are 
^ combined with the corresponding next-to-leading order results for the decay of polarized top quarks into 
■ - - dilepton, lepton plus jets, and all jets final states. As an application we predict double differential angular 
distributions which are due to the QCD-induced top quark spin correlations in the intermediate state. In 
addition to the analytic results, we give numerical results in terms of fit functions that can easily be used 
in an experimental analysis. 
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I. Introduction 



A large number of top quarks will be produced at the Fermilab Tevatron and at the CERN Large Hadron 
Collider (LHC). This makes the exploration of the interactions of these quarks one of the main physics 
issues at these facilities. Top quark spin phenomena are expected to play an important role in these 
efforts: the spin-polarization and spin-correlation phenomena reflect in detail the interactions involved 
in top quark production and decay, thus give an opportunity for precise tests of these interactions. In 
contrast to light quarks the top quark polarization/correlation effects are not washed out by hadronization. 
This is because these quarks are extremely short-lived and thus decay weakly before hadronization can 
take place [1]. 

As far as theoretical predictions on top quark pair production are concerned, the cross sections for spin- 
averaged top quark pair production have been known for quite some time to next-to-leading order (NLO) 
in QCD [2-5]. The NLO results were refined later by resummation of soft gluon and threshold loga- 
rithms; see Refs. [6-9], and the review Ref. [10] and references therein. As to top quark spin phenomena 
at hadron colliders there exists an extensive literature on theoretical investigations within the standard 
model (SM) [1 1-19] and beyond [20-32]. Yet the precision of these analyses — in particular within the 
SM — must be increased for the data samples that will be recorded at the Tevatron and at the LHC to be 
fully explored. 

In this paper we study hadronic top quark pair production and decay at NLO QCD by taking the spin of 
the top and antitop quarks into account. In particular we present the differential cross sections to order 
a^, describing the ti + X production in a general spin configuration. Short reports on parts of this work 
were given in Refs. [33-35]. 

The theoretical description of top quark pair production in proton-proton and proton-antiproton col- 
lisions at next-to-leading order in the QCD coupling involves the following parton reactions: quark- 
antiquark annihilation 

q + q^t + t, (LI) 
q + q^t + t + g, (1.2) 
which is the dominant production mechanism at the Tevatron, gluon-gluon fusion 

g + g^t + f, (1.3) 

g + g^t + t + g, (1.4) 

which dominates at the LHC, and 

g + q^t + t + q, (1.5) 
g + q^t + i + q, (1.6) 
which gives, in general, only a tiny correction. 

As mentioned earlier the top quarks decay before they can form hadronic bound states. To construct 
realistic observables, the decays of the top and antitop quarks have to be taken into account. We consider 
the SM decays of polarized (anti)top quarks both into semileptonic and non-leptonic final states, taking 
into account the order QCD corrections [36, 37]. The main SM top decay modes at that order are: 

t^bW^MVe{g),bqq\g) (1.7) 
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where qq' = ud, cs. A complete next-to-leading order QCD analysis of top quark production and decay 
thus involves the parton reactions 

gg,qq^bb + Af, (1.8) 
gg,qq^bb + Af + g, (1.9) 

g + q{q)^bb + Af + q{q), (I.IO) 

where / = ^,£,V£. In view of the fact that the total width F/ of the top quark is much smaller than its 
mass, Ti/nit = o(l%), one may expand the amplitudes of the reactions Eqs. (ll.Sll - dl.lOl) around the poles 
of the unstable top and antitop quarks, which corresponds to an expansion in powers of Tt/mt. Only the 
leading term of this expansion, i.e. the residue of the double poles, is considered here. In this framework 
the radiative corrections to Eq. dl.SI) can be classified into so-called factorizable and non-factorizable 
corrections — and, likewise, the contributions to the squared matrix elements of Eq. (11.91 ). We compute 
the factorizable corrections in the narrow width approximation Tt/nit 0. In this approximation the 
squared matrix element p of the respective reaction is of the form 



\M p oc Tr [p7?p] = pa'a^aa',pp'Pp'p- (I-H) 

Here R denotes the density matrix that describes the production of on-shell top quark pairs in a specific 
spin configuration by one of the six reactions Eqs. (II.1I) - (II.6I) . The matrices p, p are the density matrices 
describing the decay of polarized top and antitop quarks into specific final states. The subscripts in 
Eq. (II.l II) denote the top and antitop spin indices. Both the production and decay density matrices are 
gauge invariant. The production density matrices are determined from the NLO results for the reactions 
Eq. (Ot -Eq. (Ob . The decay density matrices are derived from the results for the reactions shown in 
Eq. O). 

With these building blocks the factorizable NLO QCD corrections to top quark pair production and decay 
can be computed, keeping the full information on the spin configuration of the intermediate tt state. As 
far as applications of these results are concerned, our primary aim in this paper is the study of final-state 
angular correlations of the top quark and antiquark decay products, which reflect the spin properties of 
the top and antitop quarks. For this purpose we consider the following channels: 

(1.12) 
(1.13) 
(1.14) 
(1.15) 

where I = e,/J,x, and ji, j2 denote jets originating from top and antitop decays. For these final states we 
study the following distributions at NLO in the coupling a^: 



pp.pp - 


ti+X- 


-^£+£'-+X, 


pp,pp - 


tt + X- 


-ti2+X, 


pp,pp - 


tt + X- 


-jii'-+x, 


pp,pp - 


tt + X- 





i. The double distributions 

1 da 1 



a Jcos0i<icos02 4 



1 +Bicos0i -fB2COS02-Ccos0icos02) , (1.16) 



where o denotes the cross section for the channel under consideration. Here 0i (02) describes the 
angle between the direction of flight of the lepton £+ or jet ji {i'^ or J2) in the t (f) rest frame and 
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a reference direction a (b). In particular we will discuss three specific choices for the reference 
directions a, b which allow a simple physical interpretation. 



ii. The opening angle distributions 

I do \, 

^1 - = ;r 1-Dcoscp , (1.17) 

oacoscp 2 

where cp denotes the angle between the direction of flight of the lepton (or jet ji) and of l'^ (or 
72), defined in the t or ? rest frames, respectively. 

The functional forms of the r.h.s. of Eqs. ( 11.161 ) and (11.171) hold if no kinematic cuts are applied and will 
be derived in Section IVTI 

The results presented in this work are definite predictions of QCD. Specifically, comparison of our results 
for the above distributions with future measurements will allow for detailed investigations of top quark 
production and decay dynamics. 

The paper is organized as follows. In Section |lll we present the differential cross sections for the 2-^2 
processes Eq. dl.ll) and Eq. (11.31) to order a^. In Section [nil we treat the real gluon radiation Eq. (11.21) . 
Eq. ( 11.41 ) and the qg (qg) fusion processes. In particular, we compute the differential cross sections in 
the soft and coUinear limits and perform the mass factorization. In Section |IV| we discuss the QCD- 
induced top and antitop spin effects at NLO, and we treat the issue of constructing infrared and coUinear 
safe tj spin observables at the parton level. Specifically, we compute, for all relevant parton reactions 
/ — > tlx, the expectation values of four different spin observables. These observables correspond to dif- 
ferent choices of the directions a, b used to define the angles in Eq. (11.161) . We also show that the spin 
observables have a simple interpretation as double spin asymmetries with respect to a given quantization 
axis. The spin observables serve as one building block for the calculation of the distributions Eq. ( II.16I ). 
Eq. (11.171) . The other building blocks are given in Section|Vl these are the one-particle inclusive angular 
distributions of the semileptonic [36] and non-leptonic [37] decays of polarized top quarks and anti- 
quarks to order a,. In Section |Vl| we show that the functional form of the distributions ( 11.161 ). (11.171) is 
indeed as given in ( 11.161 ). ( 11.171 ). In addition we give NLO formulae for the coefficients C and D defined 
in Eq. (11.161) and Eq. (11.171) . We further investigate the important question of whether the distributions 
Eqs. (11.161) . (11.171) are affected by non-factorizable QCD corrections. In particular we show, using the 
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results of Refs. [38-42], that the non-factorizable corrections at order do not contribute to Eqs. (11.16 



(11.171) . In Section En] we present our NLO predictions for these distributions for the dilepton, for the 
lepton-hjet, and for the jet-hjet decay channels, both at the Tevatron and the LHC. Section IVIIII contains 
our conclusions. Appendix |X| contains a collection of one-loop integrals that appear in the virtual cor- 
rections to the squared matrix elements of qq tt and gg tt, which are given in Appendices |B] and 
ICl respectively. In Appendix iDl we give the amplitudes for the processes Eqs. (II.2I) . (II.4I) . (II. 51) . and (II. 61) . 
for arbitrary spins of the t and t and arbitrary helicities of the massless partons. Appendix El contains fit 
functions for the NLO results for the expectation values of the four spin observables computed in Section 
llYl 
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II. One-loop QCD corrections to qq annihilation and to gg fusion 



In this section we present the differential cross sections of the parton processes 



q{Pi) +(l{P2) t{kuSt)+t{k2,St) , 



(II.l) 



8{pi)+8{P2) ^t{ki,st)+t{k2,st). 
at NLO in Us. Here pi, p2, ki, and ^2 denote the parton momenta, and the vectors St, Sf, with 

Sf = Sf = —I and ki ■ St = k2- sj = 0, 



(11.2) 



(11.3) 



describe the spin of the top and antitop quarks. All quarks but the top quark are taken to be massless. 
The top quark mass is denoted by m. In the calculation of the radiative corrections, ultraviolet as well 
as soft/coUinear singularities are encountered. The singularities are regulated by using dimensional 
regularization. We keep the spin vectors in 4 dimensions. Note that no scheme dependence is introduced 
in this way, in particular no '75 problem' arises. This is because the cancellation of the UV singularities 
is independent of the external spin state and the soft/coUinear singularities are cancelled in a universal 
way. 

In the (anti)top rest frame, the spin of the (anti)top is described by a unit vector (sj). We define the t 
(f) rest frame by a rotation-free Lorentz boost from the zero momentum frame of the tt quarks (tt-ZMF). 
In this frame 



As we will see in Section Hvl choosing this frame simplifies the definition of infrared and coUinear safe 
observables at the parton level. For the 2 — > 2 processes Eqs. (III. II) and (III. 21) . the tt-ZMF coincides 
with the centre-of-mass frame of the initial partons. Note that in principle the ff-ZMF is only defined 
up to an arbitrary rotation. One can resolve this, for instance, as done in an experiment where the 
direction of one of the initial hadron beams is chosen as z-axis and one orthogonal direction as x-axis. 
After having measured the 4-momenta of the t and f in this laboratory frame, the ti-ZMF can be defined 
unambiguously by a rotation-free boost. The observables that we consider in Sections |lVl and |Vl| are 
actually the same for any choice of the tt-ZMF. 
In the tt-ZMF we have 



where E and k are the energy and the 3 -momentum of the top quark in the tt-ZMF. 

In J = 4 — 2£ dimensions the Bom-i-virtual parts of the cross sections of the processes Eqs. (III. II) . (111.21) 

are, to order a^, of the following form: 




(11.4) 




(11.5) 



(11.6) 



d&{st, Sf) = do'g + da[ 




(11.7) 
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where i = qq.gg, and s = {pi + p2)^. The factors 

= ' = (A^2_l)2(j_2)2 

arise from averaging over the colours (A^ = 3) and spins of the initial partons. The 2-particle phase-space 
measure is denoted by dr2, are the Born amplitudes and My those of the virtual corrections. As 
indicated earlier, Eq. (111.71 ) still contains soft and coUinear singularities regulated in the framework of 
conventional dimensional regularization. The production spin density matrices that enter the formula 
(II.l II) are obtained from the identity 

d&{st,ST) = ^^^,dT2^-Tr [i?'(ll + Sf -x) ® (U + Sf -x)] , (II.9) 
where x, are the Pauli matrices. 



II.l. Born matrix elements 
11. 1.1. qq initial states 

For the sake of fixing our notation, we present the squared Bom matrix element for the reaction (III. II) in 
two ways, which prove useful in the presentation of the virtual, soft, and coUinear contributions to the 
cross sections. 

The squared Born matrix element reads 



where 
Here 



IfAff |2 ^ 4^2^2^^2 _ l)^^Af + C^|^s)] , (11.10) 

= 2-p2(i_3;2) _2£. (11.11) 



P=yi-^, (n.i2) 

and y is the cosine of the scattering angle in the tt-ZMF, i.e. j = p ■ k, where p is the direction of the 
initial quark (or one of the gluons in the case of gg fusion) in that frame. The spin-dependent term can 
be decomposed as follows: 

+ Ef{pi-St){p2-St)+F^\p2-St){pi-St), (11.13) 
where 





132(1 


-/)+2£, 




4 




— 

^0 ~ 


s 




_ 

^0 ~ 


--■\ 






_ 

-^0 — 


--■\ 

s 





(11.14) 
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In view of a compact presentation of the virtual corrections do"^ (see Section iri.2.11) . it is useful to write 
down the spin-dependent term in an alternative way. Defining the two unit vectors ui,U2 which are 
orthogonal to each other: 

1 



ui 



U2 



1 -y^+'fy^ 
1 

'l — y^ -|- y2j2 



1 — k — q 



where q = (p — / a/1 — and = E /m,wQ obtain 



with 



and 



Cij = -^eodij + I uuuij - -^ij + [ unuij - -?>ij ) ^2 + (whW2j + uiiUi j)e3 



e2 = 2(32(1-/) 
£3 = 0. 



(11.15) 



(11.16) 



(11.17) 



(n.l8) 



II.1.2. gg initial state 

The squared Born matrix element is given by 



where 



^2 _ ^n^a]{N^-l) 
^ ' ~ A^(1-(3V)2 



A^^(l + |3V)-2 



■SS\2 



Af + (5, • 5f) (;?i • St) {pi ■ St) + {P2 ■ St) {P2 ■ St) 



(11.19) 



+ £f (Pl • St) [Pl ■ St) +El^\y^^y{p2 ■ St) {Pl ■ St) } 



(11.20) 



and 



+ e 



(l + py)2-4 +2£^(l-(3V 



l+2|32(l-3;2)-|34^1 + (l-/)2 

1 - 2(3^ -fp^ h + ( 1 - /)2l - £( 1 - (3 V)2 - 2£2( 1 - py ) , 
32£m2 



}2,,2\ 



4(1 + P3^) 



{-|32(l-);2)+£(l-|3y)+2£2}. 



(11.21) 



The virtual corrections do^ do not simplify when using a decomposition similar to Eq. (III. 161 ). 
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II.2. Virtual corrections 



We now give our results for the order virtual corrections to the Bom cross sections of Section UTTl 
The ultraviolet singularities are removed by using the MS prescription for the QCD coupling and the 
on-shell definition of the top mass. In the following, denotes the QCD coupling defined in the MS 
scheme of Nf = 6 flavour QCD (5 massless and one massive quark), m denotes the mass of the top quark 
defined in the on-shell scheme, and /jr is the renormalization scale. The renormalized differential cross 
sections do'y still contain single and double poles in £ = (4 — J)/2 due to soft and coUinear divergences, 
which we collectively call infrared (IR) singularities in the following. These poles are removed after 
including the contributions from soft gluon radiation and mass factorization. This will be done in the 
next section. 

The absorptive parts of the one-loop amplitudes induce small polarizations of the t and f quarks orthogo- 
nal to the 2^2 scattering plane, i.e. the d&y contain also terms proportional to Zijis]p^k and Zijis\p^k , 
which are ultraviolet- and infrared-finite at this order in the perturbative expansion. These terms were 
computed in Refs. [43,44]. They do not contribute to the observables, which we investigate in Sections 
|lV|and|Vll so that we omit them in the following. 



II.2.1. qq initial states 

The NLO virtual corrections as defined in Eq. (III.7I) can be presented as follows: 



5C^ 1+132 2 /l+pj 



1 



+ —^^qq7qqdT2 , 



(11.22) 



where 



471^/ 



-EYE 



(11.23) 



Cf = {N^ — 1)/ i'2N), Je is the Euler constant, /j is an arbitrary mass scale, and 



(11.24) 



The IR-finite part y^qq is of the form 



2^,2 /»r2 



(11.25) 



with 



iTi = ^{2G+^Co{puP2AO,0) + ^^^j^s{l-^y)CQ{-pukuOAm) 



N 



IN 



(11.26) 
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The contribution G is the finite part of the sum of the gluon self-energy diagram (computed in the 
Feynman gauge) and the counterterms: 



UN-lNf 



In 



{m^)^ 6 ( 



550(^1+^2,0,0) + ^ 



— 2m^ \— 1 

{Nf-l)Boiki +A:2,0,0) + ( 1 + ^ j Boih +k2,m,m) - -Nf 



(11.27) 



The functions Co, Bq, which appear in Eqs. (III. 261) . (III.27I) . are defined in Appendix IaI Omitting the 
contributions from the absorptive parts, the term can be decomposed as follows: 



f2 



a, 

71 



(11.28) 



The 3x3 matrix C^j has the same structure as Eq. (III. 171 ) with the coefficients being replaced by e^. 
These coefficients and Ay* are listed in Appendix iBl 



II.2.2. gg initial state 

The one-loop contributions to the differential cross section of the gluon fusion reaction (III. 21 ) are of the 
following form: 

^< = ^^88{Ce^SS + ^88]dr2, (11.29) 

where the IR-singular term is given by 



^88 = ^{'^A\^F\^+\^F\^['^B + l'c+'^c\y^-y]} 



with 



'J4 



N \nN-2{Nf-l) 1 

£2 £ 6 £ 

l2 



^2 

l+P^ 

2N^ 



Cp-Nln 



N 



The infrared-finite part reads 



0~[gn In 1 ~ 



with 



^\ = s<Co{ki,k2,m,0,m) 



:i + |3^)[2 + A^^(l-[3V)] 
2N 



(11.30) 



(11.31) 



(11.32) 



+ Co{puP2,OAO)N^l+^y) 

+ Co{-pukuOAm)^[4-N\\+^y)^]{l-^y) 

+ Co{-piM.OAm)^[4-N\l-^yf]il + ^y)y (11.33) 
The term ff{2 has the structure (here too, contributions from absorptive parts are omitted) 

^2 = —{Aff+Cf!^{SfSt)+Df\{pi-S,){pi-Sf) + {p2-St){p2-St) 

+ E^^ipi ■ St) {p2 ■ St) +E^^\y^^y{p2 ■ St) {pi ■ St) } . (11.34) 
The coefficients Ay^ , . . .,Ey^ are given in Appendix ICl 

III. Real gluon radiation and (anti)quark-gluon fusion 

In this section we consider the reactions 

q{p\)+q{p2)^t{kx,St)+t{k2,St)+g{p3), (in.l) 
g{pi)+g{p2)^t{kuSt)+t{k2,St)+g{p3), (III.2) 

q{p\)+g{p2)^t{kust)+t{k2,st)+q{p3), (III.3) 
q{P\) +g{P2) t{kuSt)+t{k2,St)+q{p3) . (III.4) 

In order to deal with the IR divergences of the corresponding lowest order QCD cross sections, we 
employ the phase-space slicing method: in the case of the reactions Eq. (IIII.ll) . Eq. (IIII.2I) . we divide 
the phase-space into four regions, namely the region where the gluon is soft, the two regions where the 
gluon is coUinear (but not soft) to one of the initial-state massless partons, and the complement of these 
three regions, where all partons are 'resolved'. The Bom cross sections of the reactions in Eqs. (IIII.3I) 
and (IIII.4I) develop coUinear divergences but no soft ones. Thus, in this case, we split the phase-space 
into the two regions where the outgoing q (q) becomes coUinear to the incoming (anti)quark or gluon, 
and the complement, which is the resolved region. For convenience we parametrize in this section the 
3-parton phase space in terms of angles and momenta defined in the c.m.s. of the initial partons. 



III.l. Soft gluon cross sections 

For the processes of Eqs. (IIII.ll) . and (IIII.2I) . we define the soft region by the requirement that the scaled 
gluon energy in the cm. frame of the initial partons, Xg = 2Eg/ \fs, be smaller than some cut parameter 
.^min, where Xmin ^ 1- That is, 'hard' and 'soft' gluon refers to the decomposition 

\=@{Xg- Xmm) + 0(-^min " ■ (III.5) 
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In the soft region we use the eikonal approximation for the matrix elements M ^^'^ and ^^-^ of Eqs. (IIII.ll) 
and (IIII.2I) . and the soft limit of the J-dimensional phase-space measure. 
For quark-antiquark annihilation, we obtain 

'2 - 1 ^ ki-k2 



m 



■ + 



m 



+ 



+ 



Pi -Pi 



Kh-P^Y {k2-P^Y 

2ki ■ p2 



N 



{ki ■p-i){k2-p7>) 



2k2-p\ 



{P1-P3){P2-P3) {h-p3){p2-p3) {k2-p3){Pl-P3 



N 



h -pi 



- ih ■ P3) (Pl ■P3)^ ih- P3) {P2 ■ P3) - 

Here \ Mg'^\^ is the J-dimensional squared Bom matrix element (III. 101) . Integrating the gluon momentum 
/>3 over the soft region we obtain 



k2-p2 



(111.6) 



""soft 



^%qdT2 / \M 



qq,8\2 

soft I 



P3 



Cf 7^ — mix) In — 



2^" 



l-p2 



N^-2 1 
+ —^S{y) + -s{-y) 



2N 



(111.7) 



where 



mm 



ni-£i 



and 



s{y)=2U2 



P(l-J)- 



-2Li2 



V i-p 



(III.8) 



(in.9) 



i-p ; i-pj 

Using again the eikonal approximation the J-dimensional soft gg fusion matrix element squared takes 
the form 



m. 



gg,g\2 
soft I 



2e 



{N^ -\)M2-Mi -Mn 



P2-ki 



+ ■ 



Pl ■k2 



+ 



{N^ -l)Mi-M2-Mi2 



P2-k2 



{P2- P3)ik\-p3) {Pl-P3){k2-P3)- 

P\-k\ 



+ ■ 



{P2-P3)ik2-P3) iPl-P3){kl-p3) 



M1+M2 



Pl ■P2 



J {PI-P3){P2-P3] 
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+ 



Mi2-(A^2-l)(Mi+M2) 



m 



+ ■ 



m 



Mi+M2 + {N^ + l)Mi2 



k\ ■k2 



{ki ■P3){k2-P3) 



where 



Ml 
M2 
M12 



and |5lfl^|2 is given in Section lll.1.21 Integrating p^ over the soft region we obtain 



where 



"°soft 



Sa 
Sb 



2s 



-gg\2 



Sb + Sc + Sc 



dY2, 



N ^ Cf N 2 Cf 



+ 



■\n{x) 



N 



[A^2(i_py)+2] 



'^-2U2il-x)-'-ln\.) 



Sc = ^[4-N^{l + ^y)^] 



s{y) 



(III. 10) 



(III.ll) 



(HI. 12) 



(111.13) 



Adding the virtual and soft contributions, we obtain the following results for the remaining singular 
contribution: 



^4' + ^^soft 



m,g 



(X. S _ 1 



= -:r^CF-[3+4\n{xn 

singular ZTl £ 



singular 



1 

2n £ 



llA^-2(%- 1; 



+ 4A^ln(xn 



daf. 



(III. 14) 



III.2. Collinear contributions and mass factorization 

The order a, cross sections of the reactions Eqs. (IIII.1I) - (IIII.4I) develop collinear singularities when the 
momentum of the outgoing massless parton becomes parallel to the momentum of the incoming massless 
quark or gluon. We use, for convenience, the parameter .JCmin for characterizing the collinear regions in 
phase-space, too. For Eqs. (IIII.ll) and (IIII.2I) . we define the two collinear regions by {cosG > (1 — ;cmin) 
and Xg > Xmin} and {cos < (— 1 +Xmm) and Xg > Xmm}, where is the angle between the gluon and one 
of the initial partons in the cm. frame of the initial partons. For Eqs. (IIII.3I) and (IIII.4I) the two collinear 
regions are defined by cos0 > (1 — JCmin) and cos0 < (— 1 +-v^min)- Here denotes the angle between the 
outgoing massless (anti)quark and one of the initial partons. 
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In these regions we use the coUinear approximations for both the squared matrix element of the respective 
reaction and the phase-space measure in d dimensions. We parametrize = —z)pi + p± + 0{p^), 
p' = zpi — Pa_ + O {p\), P3- P±= 0, where Pi = pi,P2 is the momentum of one of the initial partons, p-i 
is the momentum of the unobserved coUinear outgoing parton, and p' is the momentum of parton / after 
coUinear emission. The momentum fraction z varies in the interval z E [4m^/f, 1 — 5], where 5 = x^i^ 
for qq and gg initial states, and 5 = for the gq{q) initial state. Integrating over the angles of the 
parton with momentum ps in the coUinear region, using the angular measure in J — 1 spatial dimensions, 
d^d-i = (1 -cos2 0)-^J cos0jn^_9, we obtain for the coUinear contribution to the differential cross 
sections: 



a,CF ll+z^-e(l-z)2r ^ 
^^Fe- j— ^do];'{zpuP2) +dol^{puzp2) \dz, 



dol'J{z) = -^Fe^\^j^^ + -^-z^+z-2^^^daf{zpup2)+daf{puzp2)}dz, 



(III. 15) 
(III. 16) 



doZiiz)=doZ,{z 



1 f U^ + (l-z) 
2n ^£l2 l-£ 



daf{pi,zp2) 



, ^ l + (l-z)2-ez2 
+ Cf ^ d(5f{zpuP2) \dz, 



where 



2-, 



(III. 17) 



(III. 18) 



r(i-e) 

The Born cross sections in Eqs. (IIII.15I) - (IIII.17I) are those given in Section UlTI The z-dependent terms 
in front of the JOg, which appear in these equations are, as expected, proportional to the J-dimensional 
Altarelli-Parisi splitting functions Pqq{z), Pggiz), Pgqiz), and P^g{l -z). In Eqs. (IIII.15I) and (IIII.16II . the 
integration over the momentum fraction z can be performed using the relation 



l-Xn 



dz- 



4m2 



1-z 



-£ 



1+e 

s 

In(l-z) 
1-z 



dzg{z) 



1 



1 



(1-z) 



-5(1 -z) ln(xmi„) 



+ ^e5(l-z)ln2(xmin)-fO(£2) 



The plus prescription defines distributions via 

[F(z)]+ = lim (0(l-z-ri)F(z)-5(l-z-ri) [' F{y)dy 

so that if g{z) is well-behaved at z = 1, then 

'Xz) f'g{z)-g{l) 



dzg{z) 



1 

(1-z)- 
In(l-z) 



1-z 



+ g(l)ln(l-x) 



1-z 



ri {g{z)-g{\))H\-z) ^ g(l)^^2(i 
Jx 1-Z 2 



-X). 



(HI. 19) 
(III.20) 

(III.21) 
(III.22) 
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Adding, in the case of i = qq, gg, the soft and collinear contributions, the singularities proportional to 
In(xniin) cancel for infrared- and collinear- safe observables. The singularities, which remain in do'y + 
^^soft ^^'^ ^^'coW removed by absorbing them into the unphysical bare parton density functions, i.e. 
by renormalizing these functions at the factorization scale /jf. At the level of the differential cross sec- 
tions, this amounts to adding counterterms d(5[. being composed of mass factorization counterfunctions. 
Finite Born cross sections for gq and gq scattering are obtained in an analogous fashion. Performing the 
mass factorization in the MS scheme, the counterterms read 



dafiz) 

dof{z) 
do^iz) 



1 



r(i-e) 
1 

1 



^4 



2\ e 



2\ e 



t4 



1 

£ 2n 
1 

£ 2n 

la, 
£ 2n 



Pqq {z) dof {zpuP2)+Pqq{z) dof {p\,ZP2 



Pgg{z)dof{zpuP2)+Pggiz)dof{puzp2) 



Pgq{z)dof{zpuP2)+ Pqg {z)do''^ {PUZP2 



r(i-£) 

and do^'^{z) = d(5l^{z). The evolution kernels are given by [45-49] 

z 



dz, 

dz, 
dz. 



Pmiz) 



JlN-2(Nf-\) 



(1-z) 



+ ^+z{\-z) 



Pqqiz)=CF[l^l-z) + -l'^'' 



z 



(ni.23) 

(III.24) 
(ni.25) 



(III.26) 

(III.27) 

(III.28) 
(III.29) 



III.3. Finite contributions 

In the resolved regions the four differential cross sections do'j.^^, respectively the corresponding produc- 
tion density matrices R\.^^ for the final states tt + g, tt + q, tt + q, are obtained from the helicity amplitudes 
given in Appendix O according to formula dP.ll) . In the resolved regions no singularities arise and it is 
thus possible to work in J = 4 dimensions. 



IV. Parton level results for the final states ttX 

Based on the results given in the previous sections we consider in this the inclusive parton reactions 

i^ttX, i = qq, gg, gq, gq (IV. 1 ) 

at NLO. From the different contributions JOg, dOy , do'^^^^, given in Sections|niand|lIll we can extract the 
corresponding density matrices using Eq. (III.9I) . For do'^^^^ and d&^ an analogous formula holds. These 
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and the above-mentioned matrices R'^^^ have the following structure in the spin spaces of the top and 
antitop quarks: 

where are the Pauli matrices. The functions A' determine the spin-averaged production cross section. 
The functions and B'^^ are associated with a polarization of the top quarks and antiquarks. As men- 
tioned earlier, parity invariance of QCD only allows top and antitop polarizations that are induced by 
absorptive parts of the scattering amplitudes. To NLO QCD, the corresponding polarization is orthog- 
onal to the 2 — > 2 scattering plane, i.e. the BJ and have no components in the scattering plane. The 
functions C^^ encode the top-antitop spin-spin correlations. For a general decomposition of the 2 — > 2 
density matrices and their classification with respect to discrete symmetries, see Ref. [25]. 
In the following we consider a general observable O that may depend on the spins and momenta of the 



top and antitop quarks and is thus sensitive to the different structures shown in Eq. (IIV.2I) . In particular 
we study observables of the following form 

0=4(Sra)(Sr-b), (IV.3) 

which are sensitive to the correlation C'^j^. The top and antitop spin operators are given by 

S, = ^(x®ll) (IV.4) 



and 



2 



S,- = ^(ll®x). (IV.5) 



The unit vectors a and b are arbitrary reference directions. The factor 4 in Eq. (IIV.3I) is introduced to 
have a simple relation between the expectation values of the observables defined in Eq. (HV.3I) and double 
spin asymmetries. At the parton level the following relation holds: 

The arrows on the right-hand side refer to the spin state of the top and antitop quarks with respect to the 
quantization axes a and b. A prescription of how to construct the correlated top and antitop rest frames 
in a unique way has to be given. At that point it is important to make sure that the specific prescription is 
soft and coUinear-safe. To illustrate the problem, consider the top-antitop helicity correlation, in which 
case a and b are the t and f directions of flight. An obvious frame in which these directions may be 
defined is the cm. frame of the initial partons. In fact, this is how the helicity correlation at parton 
level is usually defined in Born level calculations for hadron colliders. However, this frame can only be 
constructed by a measurement of the 4-momenta of all final state particles/jets. In particular, at NLO 
QCD, one needs to know, apart from the top and antitop momenta, also the momentum of the hard gluon 
emitted in real radiation. Obviously, this information cannot be obtained if the gluon is collinear to one 
of the initial partons. Thus the above definition of the helicity correlation is not collinear- safe and cannot 
be applied beyond the leading order. This applies also to other observables involving the momenta of the 
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top and antitop defined in the cm. frame of the initial-state partons. In the situation at hand, a suitable 
frame is the tt-ZMF defined in Section|lll As discussed in Section|ll|the top (anti)quark rest frame is then 
defined through a rotation-free Lorentz boost. We study the following spin bases, which are relevant to 
applications to the Tevatron and the LHC (see Sections IVTI and IVTlt : 

a = — b = k, (helicity basis), (IV.7) 
a = b = p, (beam basis), (IV.8) 
a = b = d, (off-diagonal basis), (IV.9) 

where k denotes the direction of flight of the top quark in the tt-ZMF and p is the direction of flight of 
one of the colliding hadrons in that frame. The direction of the hadron beam can be identified to a very 
good approximation with the direction of flight of one of the initial partons. Thus, at Bom level, the 
beam basis a = b = p coincides with the direction of flight p* of one of the initial partons in the cm. 
frame of the initial partons, and this vector is equal to the unit vector of one of the hadron beams in the 
laboratory frame. The term 'off-diagonal basis' refers to axes with respect to which the spins of tops and 
antitops produced by qq annihilation are 100% (anti)correlated [19] to leading order in a^. (For gg — > tt 
one can show that no spin basis with this property exists.) We use the definition: 

, -p+(l-y)(p.k)k ^^^^^^ 



l_(p.k)2(l-f) 

where y = E/m. The spin bases defined in Eqs. (IIV.7I) -( IIV.9I ) lead to spin observables (IIV.3I) that are 
invariant under spatial rotations of the tt-ZMF. Thus predictions involving these observables can be 
unambiguously made without an explicit prescription on how to obtain this frame. If in the case of the 
beam axis one would use a = b = p* instead, one would still obtain a coUinear-safe observable, which 
is however not invariant under rotations of the tt-ZMF. Therefore we do not consider this choice any 
further. 

In addition we find that the observable 

0' = 4SrS,-, (IV.ll) 

which is also infrared- and coUinear-safe, is also sensitive to tt spin correlations at both the Tevatron and 
the LHC. This observable can be expressed in terms of observables of the form of Eq. ( IIV.3I ). because 

3 

0' = 4£(Sre,)(erS,-), (IV12) 
1=1 

where e/=i,2,3 forms an orthonormal basis. 

As mentioned earlier, parity invariance of QCD tells us that the t and t ensembles have no polarizations 
with respect to the reference axes of Eqs. (IIV.7l )- (IIV.9l l: 

(Srv), = (Sr-v), = 0, for v = k,p,d. (IV13) 

Equation (IIV.13I ) holds for any polar vector v. These equations justify the term spin correlations for the 
double spin asymmetries Eq. (IIV.6I ). since a simple calculation shows that 

4((S, ■ a) (S,- • b)), = corr(S, • a, S,- • b),-, (IV14) 
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where for two observables Oi, O2 the correlation is defined in the standard way by 
where 



h0k = ^J{0l)-{0k)^. (IV.16) 

Using the results presented in Sections |ll| and Unl let us now discuss the observables defined through 
Eq. (IIV.3I) . Eqs. (IIV.7I) - (IIV.9I) and Eq. (IIV.I II) . The unnormalized expectation value of the observable O 
at NLO in the QCD coupling — and, as a special case, the NLO cross sections a' for the parton reactions 
in Eq. (IIV.I I) — is obtained as follows: 

= §{/^r2Tr[(i?^+i?|.+i?:„go] 

+ j dY2dzTx[{R';l,,+R[)0] + j jr3Tr[<,0]| (IV17) 

for i = qq, gg. Here dTj, is the 3-particle phase-space measure. The contributions to the left-hand side 
of Eq. (IIV.I 7 1) from the soft and coUinear regions and from the resolved region depend individually on 
the slicing parameter Xmin- For x^i^i ^ 1 the leading dependence is, by construction, logarithmic, but in 
the sum only a residual linear dependence on Xmin remains, which is due to the approximations made in 
the soft and coUinear regions. By varying Xmin between 10^^ and 10^^, we have checked for the cross 
sections and for the expectation values of the observables given below that for Xmin < 10^^ this residual 
dependence is negligible. 
For i = gq,gq we have 

o'{0), = § {/ ^r2JzTr[(i?^'fii-Fi?i.)0] + J Jr3Tr[i?Lo]| , (IV18) 

where 4>g^ = l/[2{d - 2)N{N^ -I)]- Of course, o^^ = and {o)gq = {o)gq within QCD. The state- 
ments made below Eq. (HV.17I) apply also here. 

The unnormalized expectation values Eq. (HV.17I) and Eq. dlV.lSI) are building blocks for computing the 



distributions (11.161) and (11.171) (see Section IVIl) . In computing the cross sections and these expectation 
values we have performed the phase-space integration of the resolved parts numerically. If one identifies 
the renormalization and mass factorization scales, and puts /jr = fif = /j, then the NLO parton cross 
sections and the above expectation values are of the form 

2 

&{sW) = ^[f \p) +4m..(//^)(p) +//^)(p)ln(A/Vm2))] , i = qq.gg.gq aV.19) 



O'{0a)i = ^[^fi(p) +4m.(^pj(p) +|pi(p)ln(A/Vm2))] , / = qq^gg^gq (IV20) 

and a = 1 , . . . , 4. The label a=l refers to the observable Eq. (IIV.I II) and a = 2, 3, 4 to the helicity, beam, 
and off-diagonal basis. The variable p is defined by 

P = -^- (IV.21) 
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For the beam and off-diagonal basis, one has to compute two contributions for the gq initial state: either 
the quark or the gluon direction of flight in the ti-ZMF can correspond to the direction used to define the 
axes a, b, and the contribution from hard gluon emission is different for these two cases. We have 



^(0) 

r(o) JO) 



g'gq,a — . 



(1V.22) 



The lowest order scaling functions f > and [i = qq.gg) can be computed analytically. In order to 
write these functions in a compact form it is advantageous to introduce the following functions: 

^[ln(x) + 2P], 



1 

w 

1 

w 

1 

PL 

where p = and x = (1 - P)/(l + 13). 

We then get: 



h 

^1 = o 



ln(;c)+2p + -p3 



ln(x)+2|3 + ^|33 + ^p5 



arctan ( ^— \ — (3 



(IV.23) 



/i?(p) 



JO) 



(0)( 
S^^,l(P) 

Cm 

^(P) 

4"!,(p) 



&(P) 



Tlpp 

27 
192 



(2 + p), 

[4 + p + 2p2-(16+16p + p2)£i] 



(-2 + P), 



/i^^P) 



27 

7t|3p8 + 3p + 4y^ 

~TI 5 ' 

^(0) 

Tlpp 

192 
Tlpp 

192 
Tlpp 

192 



[-28 + 5p + 2p2 - (16 + ISp + p2)£i] , 



52-7p-22p2-2p3 



+ (16 + 2p+14p2 + p3)^2 



(IV.24) 

(IV.25) 
(IV.26) 
(IV.27) 

(IV.28) 
(IV.29) 
(IV.30) 

(IV.31) 



96 248 



^ 44 



482 2 68 ./2 284 
-— p ' H 



3 , 2^4 



15 



I5P+5P 



- (16 + 34p - 64p3/2 _ _ 34p3/2 ^ j ^ p4^)^^ 



5/2. 



(IV.32) 



TtPp 



1 



28 - 35p + 14p3/2 ^30p^ - 18p^/^ +7p" +2p 



192(l + p)2 

(16 + 50p + 68p^ + 19p3 + p^)£i + 2p3/2(7 - 9p)ri 



2 10^5/2 , ^^3 , -,^4 



(IV.33) 
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The functions fj'^'^ and g^^^ that determine the scale dependence are obtained from the next formula, 
which follows from the renormalization group equations: 



7 I Pu{z) 



and likewise for the functions (p). Here, 



1 



^0 = -{nN-2Nf), 

and the evolution kernels P^j are given in Eqs. (IIIL26I) - (IIII.29I) . 

For the qq initial state we obtain the following results for the /^-dependence: 



Jqq 


Pp f 

64871 1 


Sqq,l 


_ fW 

Jqq 7 




PP f 


Sqq,2 


64871 1 




PP 


6qq,3 


324071 



Ki(2 + p)-16 



ln(x) 

T 



Ki(-2 + p) + 16|3ln(x)-y|, 

/o . . ^^ 8(-8 + p)ln(jc) 32(2p-r 
Ki(8 + 3p+4v/p) + ^ ^ ' +■ ^ ^ ' 



P 



P 



arctan 



VP 



32 1 6 
64v/pln(p) + — Vp + — 



(1) 



qq 



where 



Ki = 161n 



4p2 



127 

+ — -2A^^. 



For the gg initial state, we only give analytic results for the simple cases. With 



K2 



K3 = -31n(x)2 + 27l2-12Li2(x) 



ln(x)ln(^)-2Li2 



+ 2Li2 



'Am 



we have 



hi 



7687C 



61n 



4(3^ 



(IV.34) 



(IV.35) 



(28 + 31p)p + 3K2(32-16p + p2)+K3(16 + 16p + p2) 



^ (724 - 3328p + 7449p2) - ( 198p + 59p2 - 288) ln(x) | - 
^ ^ 61nfiV60 + 31p)|3 + 3K2(-2 + p)(-16 + p)+K3(16+18p + p2) 
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(0) 



76871 



4p2 



(IV.36) 
(IV.37) 



(IV.38) 
(IV.39) 

(IV.40) 



(IV.41) 



(IV.42) 



'88,2 



15p 
P [6 



176 - 61 18p + 7179p2) - ( 108p + 77p2 - 672) ln(x) } - 



127l2 



(IV.43) 



76871 1 (3 



In 



4(32 

2 



(60-25p + 31p2)-6K2(16-9p + 16p^ 



( 16 + 2p + 14p^ + p ~ U^^'^^^^ + 1 127p + 3399p2) 



-231p + 329p2 + 870)ln(jc) 



,(1) 



396(p2-p + i; 
p|32 



[Li2(-|3)-Li2(|3)] 



(0) 
127l2 



.(IV.44) 



For the remaining two functions gl'^JriA we give simple fits in Appendix |E| 
For the gq initial state, we obtain: 



1 



n 



87c2 192 I 3 
8p 



/-p(2-p)K2-^p(14p2 + 27p 



136)ln(x) 



1319p2-3468p + 724)|, 



135 
1 n (A 

87^2 192 1 3 

135 
1 % r4 



p(16_9p)K2--p(14p^ 



+ 



27p-328)ln(jc) 

1319p2-6258p-176)|, 

16 + 9p)K2 - |p(3p2 - 26p + 263) ln(;c) 
8p 



87^2 192 U""' " ' ''''"^ 9 
88(p2-2p + 2)[Li2(|3)-Li2(-p)] 



135 



(1439p^+1347p + 7384) 



(IV.45) 



(IV.46) 



(IV.47) 



Again we give simple fits in Appendix |E| for the more complicated functions f ^^''3 4. 

The scaling functions fj-^^ and are obtained by numerical integration. As a check we computed the 

spin-averaged NLO tt production cross sections o' as functions of p. We found excellent agreement with 

the results of Refs. [2, 4, 5]. Our result for fgq^ also agrees with the one of Ref. [2]. In the comparison 

of fj^^ one has to take into account that Ref. [2] uses a renormalization scheme where heavy and light 
quark loops are subtracted differently according to the scheme of Ref. [50], while we use the standard 

MS definition of of Nf = 6 flavour QCD. This implies that our scaling functions fj^^ are related to 
the functions (/j^'')N of Nason, Dawson and Ellis by 



/]" = (J) 



IN- 



127^2 



f- 
■'J 



(0) 



J = 11:88 ■ 



(IV.48) 



Performing the conversion we again find perfect agreement with Ref. [2]. 

For the functions fj^^ and ^-^ we give fits to the results of our numerical integrations in Appendix |El 
To account for the richer structure of the functions gj^a — compared with the results for the unpolarized 
cross sections — we have used a more general ansatz than in Ref. [2]. All scaling functions for the initial 
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Figure 1: Left: Scaling functions fqq\r\) (dotted), fql\T\) (full), and fqq\r\) (dashed). Right: The same 
for the process gg ti{g). 




states qq and gg are plotted in Figs.[T]-|51as a function of 



ri = - - 1 . (IV.49) 
P 
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Figure 3: Left: Scaling functions g"^^ 2^ (dotted), g^^^2(^) (iM^l)' S^qqii^) (dashed). Right: The 
same for the process gg — > tt{g). 
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Figure 4: Left: Scaling functions g^^^^i'^) (dotted), ^^^3(11) (full), ^f^d g'^^'^^i'^) (dashed). Right: The 
same for the process gg tt{g). 
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For large values of r|, the results for the functions g\^ and g}^ differ from the corresponding results in 
Refs. [33, 34], since in those papers the spin axes were defined using the initial-state parton c.m.s. rather 
than the ti-ZM¥. 
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V. Decays of polarized top quarks 

In the leading pole approximation for the intermediate top and antitop quarks the NLO QCD analysis 
of the reactions Eq. (II.8l) -Eq. (ILIOI) involves also the matrix elements to order of the main SM decay 
modes of the (anti)top quark in a given spin state, that is, the semileptonic modes 

t^bty^^bty^g (V.l) 

(£ = e,^, x), and the non-leptonic decays 

t ^bqq'.bqq'g, (V.2) 

where qq' = ud, cs for the dominant channels. Specifically, for computing the angular distributions (11.161) 
and (II.17I) . we need the matrix elements for the parton reactions 

i^ai+a2+X, (V.3) 

where a\,a2 denotes a lepton or a jet. Thus we require, in addition to the above production density 
matrices, the one-particle inclusive angular distributions dT /dco^Q for the decays 

t{st) a\{q\)+Xi, 

f(5,-)->a2(^2)+^2, (V.4) 

at NLO in a^. Here q\ and q2 are the momenta of a\ and a2, respectively, which we define in the rest 
frame of the (anti)top quark, and is the angle between the polarization vector of the (anti)top quark and 
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the direction of flight of ai{a2). For a fully polarized ensemble of top quarks these distributions are of 
the form 

-^ = ^(1 + k(i)cos0), (V.5) 
a cos 2 

where F^^^ is the partial width of the respective decay channel and K*^^) is the top-spin analysing power of 
ai . For the case of the standard {V —A) charged current interactions these distributions were computed to 
order for the semileptonic and non-leptonic channels in Refs. [36] and [37], respectively. ^From these 
results the corresponding unnormalized one-particle inclusive decay density matrices p and p integrated 
over the energies can be extracted. These matrices have the form 

p'ora'=^(ll + l^^'^^-qi)a'a, (V.6) 

r(2) 

pSr' = ^(ll-K^'^^-q2)a'a, (V.7) 

where qi 2 are the directions of flight of ai and a2 in the rest frame of the top and antitop quarks, 
respectively. As we work to lowest order in the electroweak couplings, F^^) = F^^^ and K2 = Ki to all 
orders in a^, if the channel ^2 +X2 is the charge-conjugate of ai +Xi. 

For semileptonic top decays, the charged lepton is the most efficient analyser of the spin of the top quark. 
In the case of non-leptonic decays t bqq'{g) this role is taken by the weak isospin % = — 1/2 quark 
from W decay. Tagging of the flavor of q or q' , besides those of the b quark, is feasible, as far as the 
dominant channels are concerned, only for q = c, albeit with low efficiency. The non-b jets in non- 
leptonic top decay may be used as analysers of the top spin in the following, more efficient way: angular 
momentum conservation and the {V — A) structure enforce that the fermion with V = —1/2 from W 
decay is, on average, less energetic than its partner with % = +1/2. This implies that the least-energetic 
light quark jet is a good top spin analyser. As far as the W boson is concerned notice that its analysing 
power is practically the same as that of the b quark: 

Kiy = -Kfo + 0(a,). (V.8) 

In Ref. [37] the coefficients K^-^) were given to NLO accuracy. In Section |VI| we need, however, the 
unnormalized density matrices (IV. 61) and (IV.7I) : that is, we require, apart from the partial widths 

r(/) = al(^ + 47ra,a5-^^ , (V.9) 

the dimensionful coefficients 

r(/) k(-^) = b''/^ + 4Tiasb[^^ . (V. 1 0) 

For the determination of these coefficients we use the Fermi constant Gf = 1.16639 x 10^ GeV^, 
m = 175 GeV, mw = 80.41 GeV, Fvk = 2.06 GeV, mt, = 5 GeV, and all other quark and lepton masses 
are put to zero. (We do not use the narrow width approximation for the intermediate W boson.) We 
obtain, for the hadronic (h) and semileptonic (si) decays t bqq' , b£v, putting the CKM matrix elements 
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\Vtb\ = \v, 



qq' 



al = 0.50645 GeV, 



a: 



si _ % 





a\ = -0.01900(8) GeV, 

af = -0.01061(2) GeV. (V.ll) 



For the LO coefficients we obtain: 



bi = 4 = ml 

bf" = -0.20663 GeV = N bQ''\ 
b^Q = -0.15819 GeV = A^Z7^ 



0' 



= 0.25840 GeV, 



bl = -0.16040 GeV. (VI 2) 

Here j stands for least energetic non-b jet, and T for an oriented thrust axis (for details, compare 
Ref. [37]). In Table[T]we give the corresponding NLO coefficients b\. The numbers in the first columns 
are for 'bare' quarks, while in the second and third columns the E-algorithm and Durham algorithm 
were used, respectively, as jet clustering schemes [37]. We demand that the four final-state partons in 
t — > bqq'g be always clustered into three jets, i.e. no jet resolution parameter is involved. This is possible 
because the leading order matrix elements are free from soft and coUinear singularities. 



Table 1: NLO coefficients b\. Numbers are in units of GeV. 





Partons 


Jets, E-algo. 


Jets, D-algo. 


Semileptonic decay 


bi [GeV] 
b{ [GeV] 
Z.?'** orb\''' [GeV] 
bl or b\ [GeV] 
b{ [GeV] 
bj [GeV] 


-0.03154(22) 

0.01354(23) 
0.00431(21) 

0.00915(22) 


-0.04218(22) 

0.01420(23) 
0.00935(22) 
-0.02336(23) 


-0.04410(22) 

0.01412(23) 
0.00852(22) 
-0.02342(23) 


-0.01074(7) 
0.00626(7) 
0.00120(7) 



Higher-dimensional distributions for t and/or t decays such as dT/ (dxid cos 0) , where Xi the scaled lepton 
or jet energy are also known. Together with the production density matrices of Sections Ull - IIVl they can 
be used for predictions of higher-dimensional distributions for the 2-particle inclusive reactions (II.12II - 

(imi . 
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VI. Angular correlations in hadronic collisions 



Let us now discuss the hadronic reactions Eqs. (II.12I) - (II.15I) . In particular we derive the distributions of 
Eq. (11.161) and Eq. (11.171) at NLO in a^, using the narrow-width approximation for the intermediate top 
and antitop quarks and taking into account only the factorizable QCD corrections. Recall that in this 
approximation the squared matrix elements of the parton reactions (II.8I) - (II.10I) are of the form shown 
in Eq. dl.l II ). We postpone the discussion of the non-factorizable QCD corrections until the end of this 
section. There we show that, at NLO in a^, they are irrelevant for the distributions defined in Eqs. (11.161) 

and (lUTb . 

With the ingredients of Sections Hn - TVl we can first determine the two-fold differential (MS-subtracted) 
cross sections 

do{i^aia2+X) 
(icos0i(icos02 

for the parton reactions / aiUi+X, where 0i (02) refers to the angle between the direction of flight qi 
(q2) of particle ai (^2) in the t (f) rest frame and an arbitrary direction a (b) that may be used as z axis in 
the (anti)top rest frame. Recall that we define the t{t) rest frame by performing a rotation-free Lorentz 
boost from the tt-7MF. For the reactions (IL12I) - (IL15I) . i.e. 

hi{Pi)h2{P2)-^tt+X ^aiai+X, (VL2) 

where h\,h2 = p ox p, the differential cross sections are obtained in the usual way: 

d<5{hi{Pi)h2{P2) ^ aia2 + X) 



JcOS0iJcOS02 



/ dXldX2 fa^{xulAf)fj^^{x2,lUF) 
„ uJ 

do{a{xiPi)b{x2P2) ^aia2+X] 
Jcos0i<icos02 



(VL3) 



where fll{x,/jf) denotes the parton distribution function (PDF) of parton a in hadron h. 

The following analysis is based on integrating over the full phase-space of the particles in the final state. 

What then is the structure of the normalized distributions 



1 d(5 

a(icos0i(icos02 



(VL4) 



if we choose as reference axes a, b one of the sets Eq. (HV.7l) -Eq. ( IIV.9I) ? Each of the contributions to 
do{i — > aia2 +X) is of the form Tr [p^^p] {s = 5, V, soft, col, c, res, depending on z), where the R'^ have 
the structure as given in Eq. ( IIV.2I ) and p and p are given by Eqs. (IV.6I ) and ( IV.7I ). respectively. All 
contributions are bilinear in qi and q2. Thus the differential cross section Eq. (IVL3I ) is bilinear in cos0i 
and cos 02. Parity invariance of QCD dictates that the top and antitop quarks have no polarization along 
the above reference axes a and b, see Eq. (IIV.13I) . This implies that in QCD 

(cos 01 ) = (cos 02) =0. (VL5) 

Here the expectation value refers to reaction (IVL2I) . i.e. (o) = J doO. (Equation (IVI.5I) applies also 
to cuts that are parity-invariant.) Thus for these reference axes the coefficients Bi and B2 in Eq. (IIV.2I ) 
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are absent in the double distribution Eq. (11.161) . and we obtain 



1 ^ a = -(l-Qcoseicose2) . (VI.6) 

o ^^coseii^cos02 4 

with / = hel, beam, off for the helicity, beam, and off-diagonal bases. Non-zero Bi 2 are generated by 
parity-violating contributions to i ttX. Computations indicate [51-53] that they are small within the 
standard model. 

Next we consider the distribution of the angle cp between the directions of flight qi,2 of particle/jet ai 
and a2, defined in the t and t rest frames, respectively. Using CP invariance of QCD and arguments 
analogous to those above, we find that this distribution is of the form 

I do{hih2^ aia2+X) 1 

^ ; = :7 ( 1 — D cos (D) . (VI.7) 

o 6?cos(p 2^ 

Defining, for ijr= /jp = m, 

Nr = -^-jY dxidX2fa'{Xi,^F)fh^{x2,IJF) 

X {,»^' +4.a.[,ii^X^)4^) , (VI.8) 



a] 1 



=2 / dxidX2 fa' {xi,pF)fh^{x2,^iF) 
a.b-' 

we obtain to NLO in the QCD coupling: 

A^i N2 N3 N4 

D= , Chel = , Cbeam = , Coff = . (VI.IO) 

a a a o 

Here //*^'^'*(p) and gf'^^\p) are the scaling functions defined in Eqs. (IIV.19I ) and (II V. 201 ). and collected 

in Appendix lEl with p = Ain^/s, where the partonic cm. energy is given in terms of the hadronic 
cm. energy y/s through 

S=XlX2S. (VI. 11) 

The total top quark width is denoted by Tj. 

Finally, we discuss the issue of non-factorizable corrections. At NLO in a^, such corrections are present 
in the one-loop contributions to the qq and gg fusion processes Eq. (11.81) and in the squared matrix 
elements for the corresponding processes Eq. (11.91) with real gluon radiation. These non-factorizable 
corrections correspond, at this order in the QCD coupling, to gluon exchange that interconnects the 
different stages of the off-shell tt production and decay process. These corrections were studied in 
the semi-soft gluon approximation in Ref. [42]. This approximation consists in taking into account 
only virtual, respectively real gluons with energies Eg of the order of Tj. This is adequate because the 
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contribution of a hard gluon to this correction is suppressed by Tt/Eg. The contributions <io^^ and do^^^ 
to the respective differential cross sections are important for calculating a number of distributions, such 
as invariant mass distributions of the ? quarks and of ?f pairs. However, the non-factorizable corrections 
do not contribute to observables, which are inclusive in both the top and antitop invariant masses [38-42]. 
A well-known example for such a case is the total cross section. 

It is straightforward to show that this result applies also to the distributions shown in Eqs. ( IVI.6I ) and 
(IVI.7I) . The non-factorizable contributions do'j^'j: and doffjr (which contain both virtual and real gluon 
exchanges with momenta being integrated over) are proportional to the respective differential Bom cross 
sections, decomposed into colour singlet and octet pieces. We represent the Lorentz-invariant phase- 
space element dT^ for the 6-particle final states in Eq. (11.81) as 

dTe oc dYttdYtdYtdM^dM} , (VI. 1 2) 

where Mf (Mj) is the invariant mass of the top (anti)quark, dTff is the phase-space element of the in- 
termediate tt state, and dTf (dTt) is the 3-particle phase-space element of t (f) decay in its rest frame, 
where the z axis is identified with the reference axis a (b) . From this decomposition, it is clear that one 
has to integrate out the top and antitop invariant masses in calculating the contribution of do'^^jr to the 
distribution (IVI.6I) . But then the theorem of Refs. [38^0] applies. 

In order to show that the opening angle distribution (IVI.7I) receives no contribution from do^^j: either, we 
recall that this distribution is due to {St ■ Sj), and St ■ Sf can be decomposed as given below Eq. (IIV.I II) . 
Thus, to NLO in the QCD coupling the above angular distributions are determined solely by the Born 
contributions and the factorizable QCD corrections. 

VII. Numerical results 

In this section we present the numerical results for the angular distributions (IVI.6I) and (IVI.7I) . In the 
absence of cuts the spin correlation coefficients can be computed in terms of expectation values: 

C = -9(cos0icos02), 

D = -3(cos(p). (VII. 1) 

In the following we present the NLO results for the coefficients C and D in the absence of phase-space 
cuts, using Eq. (IVI.8I ) and Eq. (IVI.9I) . It is convenient to rewrite these formulae: 

where / = qq,gg,gq,gq. For the evaluation of Eqs. (IVII.2I) . (I VII. 3 1) we use the NLO parton distribution 
functions of Ref. [54] (CTEQ6.1M) and of Ref. [55] (MRST2003). In addition we also use the PDFs 
of Ref. [56] (GRV98), to illustrate the dependence on the gluon distribution function. We put = 
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Table 2: Results for the quantities Z- and N\J defined in Eqs. f lV7/.2D . \VII.3\> for different PDFs at the 
Tevatron and the LHC. 



Tevatron 


/ = qq 


i = 88 


2(0) 


CTEQ6.1M 


0.032162 


3.9358-10-3 


MRST2003 


0.033193 


3.9751-10-3 


GRV98 


0.033528 


5.0460-10-3 




CTEQ6.1M 


0.010721 


-1.6864-10-3 


MRST2003 


0.011064 


-1.7534-10-3 


GRV98 


0.011176 


-2.2477-10-3 




CTEQ6.1M 


-0.016695 


2.3820-10-3 


MRST2003 


-0.017240 


2.4585-10-3 


GRV98 


-0.017436 


3.1438-10-3 




CTEQ6.1M 


0.031849 


-1.2099-10-3 


MRST2003 


0.032869 


-1.2774-10-3 


GRV98 


0.033199 


-1.6457-10-3 




CTEQ6.1M 


0.032162 


-1.2959-10-3 


MRST2003 


0.033193 


-1.3606- 10 3 


GRV98 


0.033528 


-1.7497-10-3 


LHC 


i = qq 


i = 88 


zf 


CTEQ6.1M 


0.46482 


3.2048 


MRST2003 


0.49789 


3.5616 


GRV98 


0.48075 


4.1237 


<' 


CTEQ6.1M 


0.15494 


-0.92015 


MRST2003 


0.16596 


-1.0163 


GRV98 


0.16025 


-1.1620 




CTEQ6.1M 


-0.26645 


1.3931 


MRST2003 


-0.28604 


1.5401 


GRV98 


-0.27709 


1.7621 
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Table 3: Results for the quantities z\^^ and N^]} defined in Eqs. f lV7/.2D . {IVII.3\) for different PDFs at the 
Tevatron and the LHC. 



Tevatron 


i = qq 


i = g8 


i = Q8 + 8^ 




CTEQ6.1M 


4.699-10-3 


2.312-10-3 


-3.627-10-4 


MRST2003 


4.836-10-3 


2.423-10-3 


-3.926-10-4 


GRV98 


4.854-10-3 


3.113-10-3 


-4.442-10-4 




CTEQ6.1M 


1.539-10-3 


-1.070-10-3 


1.216-10-4 


MRST2003 


1.583-10-3 


-1.143-10-3 


1.289-10-4 


GRV98 


1.589-10-3 


-1.478-10-3 


1.456-10-4 




CTEQ6.1M 


-2.643-10-3 


1.458-10-3 


-1.970-10-4 


MRST2003 


-2.721-10-3 


1.552-10-3 


-2.065-10-4 


GRV98 


-2.734-10-3 


2.004-10-3 


-2.331-10-4 




CTEQ6.1M 


4.459-10-3 


-8.647-10-4 


2.066- 10 5 


MRST2003 


4.587-10-3 


-9.246-10-4 


1.995- 10 5 


GRV98 


4.602-10-3 


-1.195-10 3 


2.258-10-5 




CTEQ6.1M 


4.515-10-3 


-9.117-10-4 


3.033 - 10-5 


MRST2003 


4.645-10-3 


-9.725-10-4 


2.975-10-5 


GRV98 


4.660-10-3 


-1.256-10-3 


3.363-10-5 


LHC 


i = qq 


i = 88 


i = Q8 + 8^ 




CTEQ6.1M 


0.04197 


1.277 


0.04565 


MRST2003 


0.04415 


1.411 


0.05120 


GRV98 


0.04194 


1.625 


0.06018 




CTEQ6.1M 


0.01329 


-0.4330 


0.03077 


MRST2003 


0.01396 


-0.4760 


0.03307 


GRV98 


0.01324 


-0.5423 


0.03082 




CTEQ6M 


-0.02417 


0.5927 


-0.07905 


MRST2003 


-0.02550 


0.6513 


-0.08540 


GRV98 


-0.02419 


0.7402 


-0.08181 
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^P = mt= 175 GeV, Nf = 5 and use the following values for af^ ^: as{mt) = 0.1074 for CTEQ6.1M, 
a,(mf ) = 0. 1062 for MRST2003, and a,(mf ) = 0. 1041 for GRV98. (The values of a, are those used by 

the different groups when fitting the PDFs.) In Table |2l we present the values of Z^^^ and A^^^^^ for pp 
collisions at ^/s = 1.96 GeV (Tevatron) and pp collisions at ^/s = 14 TeV (LHC). In the latter case, we 
only show results for the D coefficient and the Chei coefficient, since the spin correlations in the beam 
and the off-diagonal basis are very small at the LHC [35]. Table |3] lists the NLO coefficients Z^^^ and 

Nj]}. We observe the following generic features: For the coefficients Ni^a the two initial states qq and 
gg always contribute with opposite signs. For the gg initial state the GRV98 PDF systematically gives a 
larger contribution, while the other two PDFs agree quite well in most of the cases. 
Adding all contributions as in Eqs. (IVII.2I) . (IVII.3I) and taking the ratios Nr/o, we obtain LO and NLO 
predictions for the spin correlation coefficients. The results are listed in Table|4lfor the Tevatron and in 
Table |5] for the LHC, for the different decay modes of the top and antitop quarks. In these two tables 
we use the CTEQ6L (LO) and CTEQ6.1M (NLO) PDFs. As spin analysers we use the charged lepton 
and/or the least energetic non-Z?-quark jet defined in the Durham jet clustering scheme. 
At the Tevatron, the helicity basis is not the best choice, since the t and i are only moderately relativistic 
in this case. The spin correlations are large both in the beam and off-diagonal basis. In fact, since the 
results are almost the same in these two bases, the beam basis is probably the best choice since it is easier 
to determine experimentally. The QCD corrections decrease the LO results at the Tevatron. The relative 
size of the corrections varies between ~ 16% (Cbeam.off) and ~ 28% (D) for the dilepton channel and 
between ~ 27% (Cbeam,off) and ~ 38% (D) in the all-jets channel. 

At the LHC, where the gg initial state gives the dominant contribution, the helicity correlation is the 
best choice from the four distributions. Here the spin correlations in the dilepton channel are of the 
order of 30%. The QCD corrections are much smaller than at the Tevatron. In the dilepton channel 
they enhance the spin correlation by ~ 2% in the helicity basis and by ~ 10% for the opening angle 
correlation coefficient D. In the all-jet channel, the QCD corrections decrease the LO results for Chei by 
~ 8% and for D by ~ 1%. 

In Tabled we compare the NLO results for the spin correlation coefficients evaluated for the 3 different 
PDFs that we used. The CTEQ6.1M and MRST2003 distributions give almost identical results, which 
is remarkable since the individual contributions from the qq and gg initial states differ in some cases 
significantly (see Tables |2l|3l). (For previous versions of these distribution functions, the agreement for 
the spin correlations was not so striking [35].) Using the GRV98 PDF, however, leads to spin correlations 
that are up to ~ 13% smaller at the Tevatron and up to ~ 4% larger at the LHC as compared to the other 
two PDFs. This shows that the spin correlations are very sensitive to the relative quark and gluon content 
of the proton. Future measurements of these correlations may help to pin down the parton distributions 
further. 

To estimate higher order corrections to our predictions, we now study the dependence on the renormal- 
ization and factorization scales. For simplicity, we keep /j = /jf = /jr and vary ^ between nit/ 10 and 
lOmj. The results for the quantities oCbeam (Tevatron) and aChei (LHC) are shown in Fig. |6lfor the 
dilepton decay channel, where we sum over i = x. The inclusion of the QCD corrections reduces 
the scale dependence significantly. The correlation coefficients C and D are ratios in which a factor 
drops out. Therefore, at LO, the coefficients only depend on fif but not on ^r. A comparison between 
the scale dependence at LO and NLO is thus meaningless for these coefficients. Table IT] therefore only 
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Table 4: LO andNLO results for the spin correlation coefficients C and D of the distributions (IV7.6D and 
rtV7.7l ) m cai'e of pp collisions at y/s — 1.96 TeVfor different tt decay modes. The PDF CTEQ6L 
(LO) and CTEQ6.1M (NLO) ofRef. [54] were used, and /jp = ijr = nij. 







Dilepton 


Lep ton-jet 


Jet-jet 


Chel 


LO 


-0.471 


-0.240 


-0.123 




NLO 


-0.352 


-0.168 


-0.080 


Cbeam 


LO 


0.928 


0.474 


0.242 




NLO 


0.777 


0.370 


0.176 


Coff 


LO 


0.937 


0.478 


0.244 




NLO 


0.782 


0.372 


0.177 


D 


LO 


0.297 


0.151 


0.0773 




NLO 


0.213 


0.101 


0.0480 



Table 5: Results for Chei and T) for pp collisions at ^/s = 14 TeV using the same PDF and parameters 
as in Tabled 







Dilepton 


Lep ton-jet 


Jet-jet 


Chel 


LO 


0.319 


0.163 


0.083 




NLO 


0.326 


0.158 


0.076 


D 


LO 


-0.217 


-0.111 


-0.0567 




NLO 


-0.237 


-0.115 


-0.0560 
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Table 6: Spin correlation coefficients at NLO for different PDFsfor the Tevatron (upper part) and the 
LHC (lower part) for dilepton final states. 



Tevatron 




CTEQ6.1M 


MRST2003 


GRV98 


Chel 


-0.352 


-0.352 


-0.313 


Cbeam 


0.777 


0.777 


0.732 


Coff 


0.782 


0.782 


0.736 


D 


0.213 


0.212 


0.185 


LHC 


Chel 


0.326 


0.327 


0.339 


D 


-0.237 


-0.237 


-0.243 




10'^ 1 10 lo"' 1 10 

|i/m, |a/m, 

Figure 6: Left: Dependence o/oCbeam on/u — /up — l^Rfor pp collisions at \fs — 1.96 TeV in the dilepton 
decay channel. The dashed line shows the LO result, the full line shows the result at NLO. Right: The 
same for aChei/o^ pp collisions at y/s = 14 TeV. 
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Table 7: Dependence of the correlation coefficients, computed with the PDF CTEQ6.1M, on ^ = ijr= jup 
at NLOfor the dilepton decay channel. 





Tevatron 

Chel Cbeam Cgff D 


LHC 

Chel D 


mt/1 
Init 


-0.326 0.735 0.740 0.191 
-0.352 0.777 0.782 0.213 
-0.370 0.804 0.810 0.227 


0.329 -0.244 
0.326 -0.237 
0.324 -0.232 



lists the NLO scale dependence. At the Tevatron, a variation of [a between mt/2 and 2mt changes the 
results at /J = mf by ~ ± (5-10)%, while at the LHC the change of Chei is less than a per cent and D 
varies by ~ ± 3%. 

Before closing this section, we summarize how an experimental measurement of the distributions (IVI.6I) 
and (IVI.7I) that matches our predictions should proceed: 

i. Reconstruct the top and antitop 4-momenta in the laboratory frame (= cm. frame of the colliding 
hadrons). 

ii. Perform a rotation-free boost from the laboratory frame to the tt-7MF. Compute a and b in that 
frame. 

iii. Perform rotation-free boosts from the tt-7MF to the top quark rest frame and the antitop quark 
rest frame. Compute the direction qi of the top quark decay product a\ in the t rest frame and the 
direction q2 of the antitop quark decay product a2 in the t rest frame. Finally, compute cosGi = 
a-qi, cos 02 =b-q2 and cos (|) = qi ■q2. 

Note that in this prescription the t and f rest frames are obtained by first boosting into the tt-7MF. If 
this step is left out, and the t and t rest frames are constructed by directly boosting from the lab frame, a 
Wigner rotation has to be taken into account. 

The results in this paper were obtained without imposing kinematic cuts. Such cuts will in general 
distort the distributions, i.e. C and D will in general depend on the angles 0i , 82 and cp, respectively. The 
expectation values on the r.h.s. of Eqs. (IVII.ll) may still be used as measures of the tt spin correlations 
in the presence of cuts. They are then no longer directly related to these differential distributions. One 
strategy, followed for instance in Ref. [62], is to correct for these distortions by Monte Carlo methods 
before extracting the spin correlation coefficient and comparing with theoretical predictions. A future 
aim will be to directly include the cuts in an NLO event generator to be constructed with the above results 
for all relevant 2^6 and 2^7 processes. 
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VIII. Conclusions 



We have determined the differential cross sections for top quark pair production in a general spin con- 
figuration by qq annihilation, gg, qg, and qg fusion to order a^. These cross sections provide, together 
with the differential rates of polarized top and antitop decays at order a^, the NLO factorizable contribu- 
tions to the reactions given in Eqs. (IL8I) to Eq. dLlOl) . We have evaluated the corresponding spin density 
matrices for on- shell intermediate top and antitop quarks. However, our results can also be employed as 
building blocks for taking the finite top quark width at NLO into account. (For studies at the Bom level, 
see for instance Refs. [63, 64].) 

As an application of the above results we have studied a number of distributions due to top-antitop spin 
correlations. The QCD-induced spin correlations are large effects, which are easily visible in distribu- 
tions of the final-state particles. Given the size of the spin correlations, they are expected to become a 
good tool for analysing in detail top quark pair production and decay dynamics. They can be studied 
at the Tevatron and — in view of the expected large tt data samples — especially at the LHC, in the 
dilepton, single lepton and all-hadronic decay channels, by measuring suitably defined double angular 
distributions. In fact, a first measurement of spin correlations in the off-diagonal basis was already per- 
formed by the DO collaboration for the dilepton decay channel [65]. While this analysis was limited by 
the very small event sample, it clearly demonstrated the possibility of the experimental study of these 
effects. For the LHC a simulation was performed [62], showing that the experimental accuracy of a 
measurement of Chei at the LHC can be expected to be better than 10%. 

On the theoretical side, we have shown in this paper that the NLO QCD corrections to the double angular 
distributions are of the order of 15 to 40% for the Tevatron, and below 10% for the LHC. A study of the 
scale dependence of our NLO predictions indicates that the residual theoretical uncertainty due to higher 
order corrections should be < 10%. Work on soft gluon and threshold resummations will further reduce 
the theoretical uncertainties. 

In particular for pp collisions at ^/s = 14 TeV, the theoretical uncertainties are quite small. By the time 
the LHC will start operating it may therefore be expected that further theoretical progress will have 
turned top quark spin correlations into a precision tool for the analysis of tt events. 
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A. One-loop integrals 

In this appendix, we collect a number of functions that appear in the virtual corrections given in Section 
III. 21 and in Appendices |Bl|Cl Only the real parts of the integrals are given. It proves useful to express the 
integrals in terms of the following variables: 

^(1±P), 



y± 



m 

S 



1) The six-dimensional box integrals are defined by 

D^{qi,q2,q3,mi,m2,m2,m4) 
1 



^2 J (/2 _ ((/ + ^1 )2 _ fnl){{l + qi+q2)^- m^) ((/ + q^+q2 + q^)2- mj) 



The following integrals appear in our results: 
Dl{p,,p2,-k2A0Am) - 
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s{yt - 
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y- 
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>■+ ^1 


y^ 






1 



1 



-Li2(j-) + -ln^(3;-) 
--U2{y+) + Un^{y+) 



Dl{-ki,pi,p2,0,m,m,m) 



in r(i+p)(i-,)r^_^^,^^^' 



fp(l-7 



+ 2yt 



Li2(— x) — ln(.x;) In 



yt_ 

y- 



+ (p-3;)Li2 1- 



yt 



DQ{pi,~ki,p2,0,0,m,m) 



+ 



4^ 



s\2 p2(l_3;2) 



2^ 



\-y, 



. u.(,-|).u.(,-i^)]}. 



D^(pi,p2,-fci,0,0,0,m) = 
^o(-^2,Pl,P2,0,m,m,m) = 

2) The 3-point integrals are defined by 



Dl{pi,P2, -^2, 0, 0, 0, m) 
^o(-^i'i^i'i^2,0,m,m,m) 



Co(^l,^2,mi,m2,m3) = t4 / 2 



)((/ + ^i)2-m2)((Z + ^l+^2)2-m2) 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



(A.5) 

(A.6) 
(A.7) 



(A.8) 
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The following integrals appear in our results (they agree with those listed in [4]): 



Co(pi,P2, 0,0,0 

Co(pi,P2, 0,0,0 

Co{-pi,ki,0,0,m 

Co(-pi,^i,0,0,m 
Co{-pi,k2,0,0,m 
Co{ki,k2,0,m,0 

Co{ki,k2,m,0,m 

Co{ki,k2,m,0,m 

Co{ki,-pi,0,m,m 
Co(k2,-pi,0,m,m 
Co{pi,p2,m,m,m) 



Q 

s 

1 

s 



2 In (^)-^ 



+ Co(/7i,/72,0,0,0), 



1 



28^ 8 V ^ 



+ Co{-pi,ki,0,0,m), 

2- 



= Co{-pi,ki,0,0,m)\ 



1 

1 
1 

yts 



y-^-y 



1 71 

2U2{-x) + -ln^{x) + — 
I o 



+ Co(A:i,^2,w,0,m), 



-ln(;c) 
£ 

-2 In(x) In( 1 - x) - 2Li2 (x) + - In^ (x) - — 



= Co{ki,-pi,0,m,m) 



3) The two-point integrals are defined by 



1 



Bo{q,mi,m2) = 



in^J {f-m\){{l + qY-ml) 8 
The following finite parts of these integrals are needed: 



C — 
= — +5o(?,wi,m2). 



5o(^i +^2,0,0 
So(A;i,0,m 

Bo{ki- pi,0,m 

Bo{k2- pi,0,m 
Bo{ki +k2,m,m 



2 + ln(^). 
2. 

2+^ln 

BQ{ki- p\,Q,m)\y^-y. 
2 + pln(x). 



4) Finally, the one-point integral is given by 

Ao(m) = ^ / — -y 



' — = — m +m . 

m^j 8 
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B. Virtual corrections to qq it 



Here we list the coefficients A^^ and e^, defined in (IIL28I) . which appear in the one-loop contributions 
to qq tt (note that the contributions proportional to the finite part of the one-point integral (= np-) and 
those of 5o(^i 5 O7 'J^) = 2 are added to those terms that are not multiplied by an n-point integral): 



71 



2 

m 

" F 

+ P'(|3'-3) 



m 



AN VI -py 



2 



Co , ;t2 , 0, m, 0) <; (23;2|32 _ 33;2 + 21323; - 213}; + 1 ) + 

2 



+ 



A^(l-|3y) A^(l+|33;) 
5o(fci-ffc2,0,0) 



4(32 



I 1 - 3/ - _ 2133; + 43^2(32^ 



(-8); + 6p - 3p^ + 3/p3) I + - PI , 0, m 



A^2_2_ 



2A^ 



2(1 -P^ 
l-py 



(1-P3^) 



1- / 

+ —Bo{k2-puO,m) 



2(1 -P^ 
1 + Pj 



1 



2p2 

4p2(l_p2 



(2^33; + 5/p2 + p2 _ 2(33, ^ 1 _ 33,2^^ ^ 5^^-^(P - + 83;) 



A^2_2 2 
+ ■ 



,A^(l-py) A^(l + Py) 



(B.l) 



fD6(pi,p2,-fc2,0,0,0,m)A^2_2 / 1 _ p: 



TT 



2A^ 



sDl,{puP2,-hA0Am) / 1-p 



7tA^ 



+ '^Co{h,k2Am,0)lN + 2^' 



1+Py 
A^2_2 



1-pJ 



+ 



[Ni\-^y) N{\+^y) 



+ 



Bo{ki+k2AO) 
2p2 



2^.. , 3pM 5o('ti-pi,0,m)p(A^2_2)([3_3;) 



(l-p^)A^-F 



A^ 



N{l-^y) 



2Bo{k2 - pi , 0, m) p (p +3O ^ 35o(^i + ^2, m, m) 



A^(l + P>;) 



2A^ 
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+ p2^^(3P'-l)-2P'(l-P') 



+ ■ 



N{\-^y) N{l + ^y) 



(B.2) 



^2 = 



fD6(pi,p2,-/:2,0,0,0,m)A^2_2 



71 



IN 



l + P^-3p3; + 3P'3;-2py + 



2\2 



8jp(l-p2) 



m 



P2 

2p2(l-p2) 



+ ^Co (^1 , ^2 , 0, m, 0) <! Ar(2/p2 - 3/ + 2P^3; - 2P3; + 3) + 



N 



+ ■ 



+ ^"^^^^pt''^'"^ {a^(3 - 2P3. - 3/ - p4 +yY - 2P2 + 4/p2) + 1(83. + 3p3 - 3/p3) } 
5o(^i - pi,0,m) {N^ - 2) (1 - Py) 2Bo(^2 - Pi Am) (1 + Py) 



N 



+ 



5o(^i+A:2,m,m)(l+2p2)(l-/) 



2A^ 



^ I (3 - 2P3; + 2p33; - 3/ - 5p2 + 5/p2)A^ + ^ (1 - P') (P - P/ + 83;) } . 



(B.3) 



4 = ^(l-y2)(l-p2)|^Co(^i,^2,0,m,0) 
5o(^i +^2,0,0) 



-(2P + 33;)Ar+ 



8P 

A/" 



+ 



+ 



4p2 



(-2P-33; + 2P^3;)iV + 



8P 
N 



Bo{ki-puO,m){N^-2)^ ^0(^2 -;?i,0,m)p ^0(^1 +^2,m,m)3; 



1 

2p2 



2A^(1-Pj) 
A^(3p23;-33;-2p) + 



2„ o«^ , P(8-Pj) ,0^3/ A^'-2 



A/" 



A^(l+Py) 
+ 2p 



4Ar 



Ar(l-Py) A^(l + Py)) 



(B.4) 
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C. Virtual corrections to gg it 



We list here the coefficients Ay,Cy,Dy and Ey^ of Eq. (III.34I) . which are part of the IR- finite terms of 
the one-loop corrections to gg tt. 



^gg ^ |3^£'^(pi,P2,-fc2,0,0,0,m)A^2^-(l-/) [ ^ pyi^-P' 2 I 2(1-13^)^ 



271 1 1-Pj (1-N' (1-pyp 

[36+12-2[3^-3|3^ 2(4-3[3^)(l-|3^) 2(1 + 132)(1 - (3^)3 

^ l^^Dlipu -h , ,.2, 0, 0, m m), ( 1 - ,2) ^ ^^^^ + ( 1 - /) Q - /) | 

7r(l-|323;2) J 
m2Co(fci,fc2,0,m,0)A^2| 613^ - 22(3^ - Sp^ + 3 



+ <;2(2+/)|3V3(l-/)|3^-3 + 

+ 



2(34 p . ^ . -V ^ 

8p4(3-2|32) 8p4(i-p2)2 



(l-|3y)2 (l-|3y)3 
m2Co(A:i, —/> 1,0, m,m) 



p2^2 



(P'j - + |3^j' - 4)A^^ - 2|3^ - 2^y 



(- 1 + + 3p6 - 7|32)A^2 - + 4p2 + 5 - p6 2(1 - p^) (^^(1 - 1^'^)N^ + P^ + P^ - 1) 

+ 1-pJ ^ (1-Pj)2 

(-3 - 5p4 + 3p6 + 17p2)iV2 _ 2|34 ^ 2|32 - 3 - p6 4A^2(2[32^(34_5)p2 



l+pj (l + p3;)2 

^ 8A^2(l_|32-)2|32 | ^ ^2(^^^p^ „^ I ^^^^ ^^^^^^^ A^2p2 I 2 



+ 



{\ + ^yf j A^2(32 
p2(l-|32)A^3^(l-|32)(l+4p2 + |3^)A^2_|36^|32_3j34_5 

1-Pj 

(1 + P2) ((1 - |32)2A^2 ^ (34 ^ 2p2 - 5) 2(1 + p2) (1 - p2)2 
(l-py)2 {l-^yf 

^ 8^2(1^(33,) t^|-(10p); + 31-9p2)A^2^2p3; + 23-9p2 

(l-p2)(|32-5)(A^2_i) (P4 _ i3p2 ^ 18)^2 _ 14 _ p4 ^ 9p2 

^ l-2P); + p2 + rr^^ 

8(l-p2)((3p2-5)A^2_3(32^3^ 8(1 -p2)2((l -(32)A^2^(32) 

+ 



(l-py)2 (l-(33;)3 
^o(fci + fc2, 0, 0)A^'( 1 - /) (-2 + (-3 + 4p2 + 2p4) /) 

4(l-p23;2) 
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+ ^4l^\|!'pT;2f {2PV(1-PW + (1 -/)(2-py-2pV)A^2 + 4(l -/)} 
+ I (-60|3 V + 23|3^ + Sepy - 41p4 - Tlp^ + 36) A^^ 



24p4(l - (4(2 - p2)iv4 ^ ((34 ^ 6(32 _ 13)^2 _ p4 _ 3p2 ^ 5) 



+ 
+ 



48p^ (A^^ - 1 ) ( 1 - P^)^ ( ( 1 - P^)A^^ + P^) 

12p^((2p^° + 17p^ + 3 + 2p^ - 1 15p^ + 43p^)iV^ + 24p^ - 4p6 + 12 - 2p^° - 6p^) 

(3 + p2)(l-P);) 



, 3pnA^^-l)(l-p^)(p^-5)(iV2(3 + p2)-4) , 4p6(l-p2)A^3^y ^ , ^ _ 
^ (3 + p2)(i + p2 + 2pj) ^ (1-PV) ^ ^ ^ 



^ IDl^jpu p2,-k2 AO Am)NW f 10-3^^ 6(1 -p2) 2(1 -P^ 
^ 71 1 1-py (1-Pj)2+ (1-Pj)3 

-p^jS + 2p33; - py + 4p2 - 4P3; - 4 ■ 



+ 
+ 



l-p2 

^§(-^l,Pl,P2,0,m,m,m)m2(2+iV2(l-p>;))(l + p2) r8-p2 6(1 -p^) 



7rp2yv2(l+p);) il-Pj (l-p3;)2 



2(l-p2)2 -p5y + p4 + pBy _ p2^2 _ _2^y-4 
+ (l-py)3+ l-p4 

4Dl{pi,-ki,p2,0,0,m,m)m^ f 2 2p4 + 5 2(l-p2)(2 + p2) 2(l-p2)2 



7t 1 1-P2 ' l-py (1-Pj)2 (1-Pj)3 

+ Co(^i,^2,O m,0)^2^2 . _ ^^^^^^^^^ ^ ^^^^ - 7)p2 + 3 - ^ + - ^^P' - 



2P4 1^ V . .r- . V . ' 1_P^ 

16P^(l-p2) 8P^(l-p2)2 
(1-Py)2 + (1-Py)3 

+ ^o(/m ^ -in,0 nu,nW f ^^3^ _ ^ ^2^2 ^ , _ p4)^2 _ ,p2 _ 



P2^2 

(l-p^)((l + 3p2)A^2_3_p2) 2(l-p2)2(2A^2p2_|32_i) 8jV2(i_[32)2[32 

I-Pj + (1-Py)2 + (1+Py)3 

(3p6+ llp2 -5p4 - l)Af2 _ p2 _ p6 _ p4 _ 1 4^2(1 _ (32^(3 ^ ^2^(32 



l+pj (I + Pj)^ 

P2]V2 
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- |3^)A^3 ^ ^[34 _ [32 ^ p6 _ ^^^2 ^ 3 ^ [32 ^ 3|34 ^ pe 
(l-|34)(Ar2(i- [32) -3-132) ^ 2(1 + |32)(1-|32)2 

+ 



(1-Py)2 (1-Py)3 
5o(^i - PI , 0, m) (2 - Ar2( 1 + py)) (1 - p2) r (_9|32 ^ ^ ^ ^ iopy)iv2 - 2Py - 19 + 9p2 



8A^(l + p3;) \ 1-p 



(l-p2)(jV2-l) 4(p2-8)(jV2-l) 24(l-p2)(jV2-l) 8(l-p2)2(jV2-l) 

+ ^l-2P3; + p2 + l-Py + (1-P3;)2 (1-P3;)3 

Bojki +k2,0, 0)N^ (1 - y^) (2 + (-3 + 4p2 + 2p4) /) 
4(l-p23;2) 

2pV(l-p2)2A^3^(-2-(2 + 5p2 + 2pV + P^(l+2p^)/)A^^ 



Bo{ki+k2,m,m) 

1 - p2 r (60pV _ 43[34 + i20p2 - 23p6 - 36p23;2 _ 3e)N^ 



l-p2 

3p2(-9p4 + 31p2 + 4p2);2 _ 4)A^2 ^ i2p4 - 48p2 4A^4^6(36 _ 26p4 _ 2ip2 + 9) 



■ + 



+ 12p 



1 - p2 1-p); 
2 (4p^ + 13p^ + 54p2 - 3 + 2p^)iV2 _ 2p8 - 12 - 38p^ - 6p^ - 12p^ 

(3 + p2)(l-P);) 



3p^(l-p^)(iV^-l)(3iV^+iV^P^-4) 24p^( 1 - p2) (N^-l) {AN^ - p^ - 3) 
+ (3 + p2)(l + p2 + 2P3;) + (1-P3;)2 

48p4(l - p2)2(Ar2 _ 1)2 A^^N^Nfj^ 



{\-^yf 1-P23;2 



+ {y^-y). (C.2) 



(l_[32)/^o(i'i'i'2,--t2,0,0,0,m)A^2p (32 



7tp2 



+ 



^§(-^i,Pi,P2,0,m,m,m)(2 + A^2^1_(3_^)) 



7tp2iV2(i^[33;) 



2p2 



+ 



2py 



l-py (l-py)2 l-3;2 



+ 
+ 



2(l-p2)2D6(pi,-fci,p2,0,0,m,m) 
7ip2(l-3;2)(l_|323;2)2 

Co(^l,^2,0,m,0)jV2 [^^^ , (3 + P^)(2p2 - 1) 4p2 

4j34 -^U P J + 

Co(A;i,-pi,0,m,m) 



4p2A^2 



l-Py \-y_ 
( 1 - p2) (2Ar2p2 + p2 _ 4) 2(1- p2) (Ar2 - 1 ) p2 



1-py 



2(2-3p2 + p4)A^2^(34^3(32^4 4(1 _ [32)^2(32 (x_^2^^2 j^2^y_2 



+ 



1 + Pj 

Co{pi,p2,m,m,m) 



+ 



-4- 



1-/ 



[32^2 



A^2 + [32 + l 1_[32 



l-pj (l-p3;)2 1-3; 
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+ 



+ 



+ 



+ 



Bo{ki-pi,0,m){2-N\l + ^y)) 



(l+A^2)(32 



+ 



Bo{ki+k2,0,0)N^ 



3-2132+ 



4|32-3" 



2p2iV2i 



-A^^(3A^^(1-P^) + |3^) + 



(1-P2)(1_P3;) -(l-p2)(l-3;) 

4(Ar4-Ar2 + 2)p2 



1+y 



(5p2 - 3)A^4 _ ^2 ^2^4 ^ 5p2 _ 1)^2 ^ 5[34 ^ [36 ^4p2 
p4(4A^4 _ 2^2p2 _ 3 2p4(l - p2) (N^-l) (2A^2 _ ly 



(i-Pj)2 



(C.3) 



^ A^^^^(Pi.P2,-^2,0,0,0,m)r p,_^py_^p^_^^_^_l+2p, + p2 7P2-19 



( 1 - P2) (9 - P2) ^ ^ (l-p2)2 I ^ iV2^^(^^^^^^_^^^o^o^o^^) r _ ^ ^p2^2 _^ ^ 



(1-P);)2 (l-py)3i 71 

^l-p2 ^ (l-p2)2 ^l+2P3; + p2 



1+pj (l + p3;)2 p2(l-j2) 
D6(-/:i,pi,p2,0,m,m,m)(2 + A^2^1-p3;))/^p2_^p^_^^(l-p2)(3p44-4p2 + 9) 



_ ^ (l-p")(l-p2) (l + p2)(3p-^ + 6p2-l) ^ ,(l-p2)P3; 



(l-p3;)2 l + Py l-3;2 

^ ^^(-^2,Pi,P2,0,m,m,m)(2 + A^2^1^(3y)) ^^^^ ^ ^ ^ ^ (3_p2)(i_(32)2 



47Cp2A^2(i^P^) I K . . (1+Py) 

(l + p2)(p6-15p4-13p2-5) p(l+p)^ (l-p)3p ) 

(l-p2)(i-P3;) ^ (l-p)(l-y) (l + p)(H-j)J 

2D6(p,,-fci,P2,0,0,m,m) f 2p4 + p2 + 9 (1 - p2)(5 + 3p2) ^(l-p2)2 



71 1 1-Py (1-Py)2 (1-P3;)3 

(l-p2)(3 + 2p2) (l-p2)2_ i+2Py+p2 



1 + Py (1 + Py)2 p2(l-3;2) 

(l-p2)Co(^i,^2,0,m,0)Ar2/ ^^^_p2^^^^p^^^l5p6-13p2-24p4 + 6 



4p4 I V K.v . K.. (I_p2)(i_py) 

+ 4P^,AZ|L_16P^^ + P2l±lf -4P^^ + 8P2 + 



(1-Py)2 - (1-P3;)3 - I+P3; - (1+P3;)2 - (l_p2)(l_3,2) 

^ (l-p2)Co(fci,-pi,0,m,m) f ^^2(1 I I p2) ^ (2 + 3p2 + 3p>2_5_2p2_p4 



2p2A^2 I V . . K / (1-pj) 
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2(1 - p^) (2 + - 3^^N^) + 3p2)^2 + 1 _ 4[32 _ p. 



4p2iV2 ( 1 _ (32) ^ ^2 ( 1 ^ 2P3; + _ 2( 1 + Py) 



(1 + Py)2 

+ ^°(^\7,^'^'-'-U(l-P^)(3 + 2p,-p>^+4(l-p^)-3-^^-P')^(^-^'P') 



2p2^2 |v K /V ' K / 'V (1 + Py) 

(5 - 3p2) (1 - p2)2jV2 + 5 + 2p2 + 9P^ 4p2jV2(l - p2) (3 + p^) 

1-Py (1-Py)2 
16p^A^2(l-p2)2 _ ( 1 - p2) ( 1 + 2Py + p2) A^2 _ 2 _ 2p(3y + 3P + P^j) 



(1-Py)3 l-3;2 
Co(pi,p2,m,m,m)(l-p2) [ 4(2+iV2)(i ^2p3; + p2) (l-p2)(3-p2) 



2A^2(32 1^ (1_(32)(1_3;2) (1+P)^)2 

p2( 1 - p2)2A^3 _ (32 ^ 1 _^ p2) (3 p2)^2 _ 7 _ 3p2 _ _ p6 ^ ^ ( ^ _ (34) 



+ 



(1-P2)(l-Pj) (l-p).)3 

5 + 3p^-2iV2(i_|34) (Ar3(32(-i _ (32) ^ ,-2 _ 3^2 _ (34)^2 _ ^ _ 5(32 _ (34-) 

(l-N^ ^ (1 + Pj) 

5o (fci - i^i , 0, m) (2 - A^2^ ^ (3^) w ^25 - 2p2 + p4)A^2 _ 17 _ (34 _ 6[32 



4Ar25 \ 1 - Py 



2(l-p2)((9 + p2)iV2-p2-5) 4(l-2jV2)(l-p2)2 32(Ar2-l) 



(1-Py)2 (1-Py)3 (3 + p2)(l-2Py+p2) 

(-25 + p4 - 9p2 + p6)jV2 + 1 - 23p2 - 9p4 - P^ jV2(l + Py)-P(P + y) ) 
(3 + p2)(l + P3;) + p2(l-3;2) / 
\^-m \, ^^2 (5A^2_3)(i_p2) ;v2(l-p2)2 



go(^2-pi,0,m)(2-Ar2(l-Py)) / _ 2 

iV2f(l_P3;) \' 



+ 2 



l + Py (1 + P3;)2 

Ar2(l+2Py + p2) + p(p2y + 2p + y) \ 
p2(l-y2) / 



+ Mi±^^/6(i + p2)(i_2p2)(l + p,)-.(^-llP^-l^P^ + 4P') 



+ p 

+ 



2p4f 1^ V ' M K 7V ' r-^y j_p^ 

22(1 - p2)2A?3 _ (3 ^ 5p2)Ar2 _ 8 8(A^2(32 + 2 + yv2) ( 1 _^ 2p); + p2) 



(l+pj) (l-p2)(l-3;2) 
(-2p2(l - p2)3A^3 ^ (2 + 1 lp2 + 24p4 - 5p6)A^2 + 8(1 + 3p2)) 

(1-P2)(1-P3^) 



+ ^|6^^(l + P.)((3pV4p2-3)^2_p2_p4)^(32P^-107P^-39p2+18)A.4 
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+ 



3|3^(39|3^ + 3|3^ + 2|3^ - 2)N^ - 3|3^(8 + 27|3^ + lOp^ + p^) 

6|34(1 - (32) ((|32 - l)N^ + (|32 + 9)A^2 _ p2 _ 3) ^2(34(1 - p2)2(^2 _ ^2A^2 _ l^ 



+ 



(l-|3y)2 (l-(33;)3 
;i - p2)(32^4(32[32 3(1 _ p2)(34(((32 ^ 5)(2(32 + 7)A^2 _ p4 _ ^4^2 _ 35) 



l + Pj (3 + |32)(l + (33;) 

6N^^\\-^'^f{N^-\) 2A{N^-N^ + 2)p2(i + p2 ^ 2|3y) 
(l+|3y)2 ^ 1-/ 

+ (3 + |32)(l-2|3j + |32) r^pv /■ ( • ) 

D. Helicity amplitudes for the real emission processes 

In this section we give explicit results in terms of helicity amplitudes for all real emission processes 
needed to compute the 2^3 spin density matrices i?5-es' ^ — I^^SS^IS^^S^ discussed in sections HiH and 
Hvl Generically, the matrices R'^^^ can be obtained from 

£|'r|2ocTr[7?i,,(ll + Srx)®(ll + Spx)] , (D.l) 

where the amplitudes T are given below and the sum stands for averaging or summing over colours and 
unobserved spins of the initial- or final-state particles. For the evaluation of the spin correlation observ- 



ables introduced in Eq. (IIV.3I) one can also use the right-hand side of Eq. (IIV.6I) . 



D.l. The process gg ttg 

In this subsection we present the Bom amplitudes for the process 

g{pMP2)g{p3) t{h)t{k2) (D.2) 

for all possible gluon helicity, and for t and f spin configurations in 4 dimensions. The gluon 4-momenta 
in (ID. 21) are incoming, and the tj momenta are outgoing. The amplitudes for the process gg — > ttg 
can be easily obtained by crossing from the amplitudes given below. It is convenient to decompose the 
amplitude for (ID. 21) with respect to its colour structure. Performing this decomposition the T matrix 
element is then of the form: 

'^{888-^tt)= {TaiTajTa^)c,CTA{pi,Pj,Pk,kl,St,k2,St), (D.3) 

{/j,^}e/'( 1,2,3) 

where the sum runs over all permutations of the indices {1,2,3}. The Ta are the generators of colour 
SU(A^) in the fundamental representation. The indices a/,Cf , and cj label the colour of the gluons, quarks, 
and antiquarks, respectively. The amplitudes A(;?i,;?2,i?35^i7>^f,^27>^f) are the so-called colour-ordered 
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subamplitudes. They can be computed in a very compact way by using spinor helicity methods [57-60]. 

2- 



We use the following notation [57] for massless spinors u,v of helicity ±i: 



w(p,±) = v{p,t) = \p±), 

u{p,±) = v{p,T) = {p±\. (D.4) 
For the spinor products we use the short-hand notation: 

(PiPj) = {Pi-\Pj+), 
[PiPj] = {Pi+\Pj-), 
{pi±\pj\pk±) = {pi±\i/j\pk±) = {pi±\pjT){PjT\Pk±), (D.5) 

for pj = pj = p^ = 0. The polarization vector of a gluon with momentum k can be written in the 
4-dimensional helicity scheme as follows: 

The momentum q denotes a 'massless' (q^ = 0) reference momentum, which is otherwise arbitrary. For 
a top quark with mass m and spin vector 5, we use the formalism of Ref. [58]: 



m 

with 



u{khSt) = —^\qj) + \q2), (D.7) 
ki = qi+q2, and St = {q\ - qi) / m (D.8) 



and 

v(^2,^f) = ^k+)-|r2), (D.9) 



m 
with 

k2 = r\-\-r2, and sj = {r2-r\) /m. (D.IO) 
Note that qi,q2,r\, and r2 are massless momenta with 

I{qi-q2) = 2{n-r2) = n?. (D.ll) 

We use the short-hand notation 

{PiM%Pj\Pk^)=%{PiMPj\Pk±)- (D-12) 

(If the momentum of the top or antitop quark is sandwiched between two spinors one must first de- 
compose it into the sum of massless momenta qi -[■ q2 or ri -|- r2.) For specific gluon helicity states we 
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obtain: 



[piq\\[p\P2][p2P3\[p3q\] [p\q\\[p2q\][p3q\]^ 

4(Pl-^l)(P3-^2) ^^^^l^^l^^^^ [(P2?2)(?l+|fcl|pi + ) - (PlP2)(?l+|pik2+) 
(^1^1/^1 +P2|?2+)(?l"^|/i2|p3^) 



2(p3-^2) [P\P2] 
1 1 



2(pi - A:!) [p2Pii\ 
A(l+,2-,3-) = - 



(?l'^|P3+P2k2+)(?l+|A:i|pi+) + (?i+|pi|?2+)(?l+|p2 + P3|pi+)]), (D.13) 

+ ry- ^-^ _ 1 (P2P3)((n^2)[Pl^2](Pl + |P2+P3k2 + ) ^ W2) [Pl^l] + |P2 + P3 Ih^)) 
2 (m2 + (fci ■/:2))(/?ip2)[piP2][P2P3]('"l'"2) 



^1 {Pl^\h\p2'^) / 2 + I P3|r^+) 

4 (pi •^i)(pi7?2)[j:?ij:?2][piP3] MP2P3] ('"1^2) 

1 (?1?2) 



(P3-^2) (n^2) 



[Pl?l]((P2P3)(P3'-l)blP3] + (P2ri)(pr|/r2|p3 )) 



1 2 

-(P2?2)[Pl^2](;?l + |^2|P3^) - 7 t[P1'"2](P1^|P2+P3|'?2^) ), (D.14) 



(P3 • ^2) [P2P3] 

2+ ^ 1 (piP3)^([P2r2](g2~|7?l +P3|P2~)(rir2) - {q\q2)[p2q\]{ri-\pi +P3\p2~)) 
2 (P1P2)(P2P3)[P2P3](^1^2)[P1P2](^1 -^2)+"^^ 

+ ^7 w u ^'^if^^ u TT f(^l^l) (^3^2) [P2?l] [P2^2] {nr2) 

2 (PlP2)(P2P3)iP2P3j('"l'"2)iPlP2j(Pl -^l) ^ 

+(pi?2) [p2^2] {p-i^ \p\ - q2\p2^){nr2) + {p-in) {qm) \p2q1] {pr\qi + ?2|p2~)) 

-(P3~|^2|P2~)((P1?2)[P2^2](^1^2) - [P2?l] (?1?2) (pi^l))) , (D.15) 

2- 3+N ^ 1 (PIP2) ([P3?'2] (g2~ \P1 +P2\P3~){rir2) - {qm) [P3qi] {n' \PI +P2\P3~)) 
2 \p2P3]{P2P3)\piP2]{rir2)iki-k2) + rrfi 

1 {P2~\k2\P3~){[p3r2\{q2~\Pl+P2\P3^){rir2) - (gl'?2) [P3'?l] (^l" |Pl +P2|P3~)) 

2 [PlPi] {P2P3) \P\P2] {r\r2) \p2P3] {P3 ■ h) 

+ {pC\h\P3^){[P3r2]{P2q2){nr2) - {P2r\){q\q2)\p3qi\))- (D.16) 

The subamplitudes for the remaining gluon helicity configurations can be obtained by exploiting, for 
example, CPT invariance. Calculating the squared amplitude in terms of colour-ordered subamplitudes, 
we have: 

t2_ 



W{ggg^tt)\^ = N'^^l £ \AiiJ,k)\' 

^ {i,j,k}eP{l,2,3) 
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{;-,;,fc}eP(l,2,3) 



{/,j,A:}e/'(l,2,3) 



>.(D.17) 



The structure of (ID. 171) must be the same as that of the T matrix element squared of the reaction Z, y — 
ttggg, which was computed in Ref. [61]. We find agreement with this result. 

D.2. The process qq ttg 

Here we give the Bom amplitudes for the reaction 

q{pi)q{p2) t{ki)t{k2)g{p3). (DAS) 

The colour decomposition of the corresponding T matrix element is given by 



T{qq^ttg) = ^5^^(r,)„-Ai + i5„-(7;)?^A2 + 5^K7:.)f^A3 + 5,^(7;)^,-A4, (D.19) 



with q,q denoting the colour indices of the massless quarks, while tj denote the colour indices of the 
top quarks. The amplitudes for the processes 

qg tiq, qg ttq (D.20) 

can be obtained by crossing from (ID. 191) . Using the above colour decomposition, the squared matrix 
element is given by 



) 



where we have used 

Ai+A2+A3+A4 = (D.22) 

in order to simplify the term that is subleading in the number of colours. Using again the 4-dimensional 
spinor helicity methods as in the previous section, we obtain the following results for specific helicity 
configurations of the massless partons: 



M{pi ,P2 ./'3+) = +77 — ^- — \x f- / ^ A P3'^\h\pi^){piq2)[p2r2] 

^[Pi-P2) [Pi,-h)^ {piP3) 

1 1 , , ^ 1 



+7 ri r{P3 \ki \pr) (pi^i) (^1^2) [P2qi] - 7 r (i^i^i) (^1^2) [P2P3] [PsQi] 

{PlP3) \r\r2) {nrj) 

+77 77 -(Pl^2)[;?2;?3][;?3''2] + 7 r{p3'^\k2\pi^){piq2)[p2r2] 

'+{Pl-P2) {P3-k2)^ {PiP3) 

^ ^ {p3^\k2\Pi+){piri){qiq2)[P2qi]), (D.23) 



{PIP3) {nr2) 

M{pr,P2^,P3^) = +77 — , 2 (7 ^7 y{PiP2){Pin){qiq2)[P2P3][P2qi] 

4(P2-P3) {ki-k2)+m^\{pip3) {nr2) 



Al 



+ (P1?2)[P2P3][P3^2] - J^^{PlP2){piq2)[p2P3][P2r2\ 



^-^^in){qm)[p2P3][pm]), (d.24) 



^3(P1 ,;?2 ,7^3 ) = +T7 TTl , N . 2 ! ~7 W rWl^2)LP2;?3 LP2^1 

4(pi-P2) (^l-^2)+m^ V {PlP3){rir2) 

- (Pm) \p2p3] [P3r2] + (P^P^) (^1^2) [P2P3] \p2r2] + -j^^ {Pin){q\q2) \p2P3] \p3q1] 

11 1/1 

+ T 



4(P1-P2) (p3-*;i) ^ {P\P3) 

11 1 

{P3^\h \pi^) {piri){qiq2) [p2qi] + 7 r{p\n) {qm) \p2p3] [P3q\] ] , (D.25) 



(pips) {nri) (rm) 

Al(pr,P2-,P3-) = +^(^(^(-(piP3)(p3^2)[p2r2] 

+r— ^(^2+|A:i|p3'*")(pi?2)[p2r2] - . . , ^ (p2"^|A:i|p3"^)(pin)(gig2)[p2gi]) 

+77 si TT ( - 7 r{PlP3){P3n){qiq2)[p2qi] 

4(P1-P2) (P3-^2) ^ (n^2) 

-|^^(P2+ |fe|p3^) (Pl?2) [P2/-2] + (;:^^^2+|^2|p3 + ) (PIH) (?1?2) [P2?l]) , (D.26) 

A2(pr,P2~,P3~) = 77 — ^- — — , I . a f' r ^ 1/ ^ \ (piP3)(pi?-i)(gig2)[piP2][p2gi] 
4 (pi -^3) {ki-k2)+m^ V [P2P3J (n^2) 

+ (PlP3>(P3?2)[P2^2] - -^^:^(PlP3)(P3n)(?l?2>[P2?l] 

+ |^^(PlP3)(pi?2)[piP2][P2r2]), (D.27) 
A3(pr,P2",P3") = 77 j\. 2 (- r ^ 1/ ^ \ (PlP3)(pin)(glg2)[piP2][p2gl] 

4 (pi • P2) (^1 • k2) + \ [p2P3j (^1^2) 

-7 T (PlP3) (P3n) (qm) \p2q\] + 7 T (piPs) (Pl?2) \P\P2] \p2r2] + (P1P3) (P3?2> \p2r2]) 

{rir2) [P2P3J ^ 

+ (P1P3)(P3?2)[P2^2] - 7 T(P2"^|fcl|P3"^)(pi?2)[p2^2] 



2(pi •P2) 2(p3-/:i) V [P2P3] 

1 1 

+ 7 Tl \(/'2+l^lb3 + )(Pin)(?l?2)[P2?l]) 

[P2P3j {nr2) ) 

+77 777 — , 2 ( ~ kP\P3)kP3q2) [P2r2] + 7 r(PlP3) (j>3n){qm) [P2^l] 

4 (pi •P3) (^1-^2) V (rir2) 

-T T (PlP3) (Pl?2) \p\p2] \p2r2] + 7 T 7 r {p\P3){pm){q\q2) [p\P2] [p2?l] ) • (D.28) 

[P2P3J [P2P3J ('"1^2) ^ 

The remaining two helicity configurations can be obtained by exploiting the parity invariance of QCD. 
The amplitudes A determined in this way may get an additional phase, which cancels, however, when 
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calculating (ID.21I) . 



E. Fits to NLO results 

In the following we give results in the fomi of fits for those functions discussed in Section|lVlthat are not 
given in analytic form. These fits can be used to obtain predictions at the hadron level by convoluting with 
the corresponding parton distribution functions. At small (3 values, the accuracy of the fits is at the per 
cent level. For larger (3 values, they are less precise — but still precise enough for all phenomenological 
applications. Note that because of cancellations between different contributions, the 8 digits of precision 
of the parameters a, are sometimes important. 

For the scale-independent part in the reaction gg ttX we choose an ansatz similar to that used in 
Ref. [2]: 



.(1) Ji) 



2^2 366 2^ , 11 2 



12pln(8p^)^-— pln(8p^) + -x 



153671 

+ (3[ai +|32(a2ln(8|3^) +^3) +a4|3^1n(8|3^) +p^(a5ln(p) +a6ln(p)^) 

+ p(a7ln(p)+a8ln(p)^)] 

+a# + aiop^ln(p)+p3(aiiln(p)+ai2ln(p)2). (E.l) 

The coefficients c^^./ are determined from the behaviour for small p. Fitting this function to our theoret- 
ical prediction, we obtain the values given in Table [51 
For the reaction qq — > ttX we chose an ansatz similar to that of Ref. [2]: 

f^ls^l = c,.v4[yPln(8P^)2-fpin(8p2)-^; 

+ |3p[ai + |32(a2ln(8|32)-fa3)-M3^(a4ln(8|32)+a5) 

+ fl6P^ln(8p2)+a7ln(p)-t-a8ln(p)2]. (E.2) 

The fitted parameters are given in Table|5| For the reaction qg ttX one has to consider two possibilities: 
either the quark or the gluon direction of flight in the tt-TMSi corresponds to the direction used to define 
the axes a,b. Li the case of the beam and off-diagonal bases, these two choices give different results at 
the parton level for the contribution from hard gluon emission. We denote the scaling functions for the 

first case by a and for the second case by g'"^}-:, . . As fit ansatz we choose 

/il^4gL,. = |3[|3'(ailn(|3) +a2) + |3^a3ln(|3) +a4) +p2(a5ln(p) +a6ln(p)2) 

+ p {aj ln(p) + as ln{pf)] + [ag + aio In(ri) +an H^f] 

+ p3[ai2ln(p)+ai3ln(p)2]. (E.3) 

The fitted values are shown in Tables flOl and fTTl 

^From the analytic expressions for the leading-order results, the scaling functions determining the /j 
dependence can be obtained by a simple convolution with the corresponding Altarelli-Parisi evolution 
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Table 8: Fit parameters as defined in Eq. \E.1\> for the functions fgg ,gggl 





Ji) 
Jgg 


(1) 
^ggA 


1 

^gg,2 


(1) 
^gg,i 


(1) 

^ggA 




1 


-3 


1 


-1 


-1 


a\ 


0.10981761 


-0.31527984 


0.10618585 


-0.10654108 


-0.10708321 


(32 


-1.0630975 


-0.27483778 


-0.1722542 


-0.091362874 


0.062373685 


a-i 


-1.2553496 


-0.98582485 


0.053571485 


-1.4480832 


-0.31938862 


a4 


-9.4872673 


-4.7958679 


0.46063565 


-1.7773284 


0.34555702 


as 


-3.6880817 


-1.9966529 


0.40949176 


-1.5264395 


0.038824512 


ae 


-5.4736051 


-2.8306013 


1.3268284 


-1.3615869 


0.052584228 


a-] 


2.0430787 


1.0567063 


-0.66333002 


0.12455129 


-0.35180567 


as 


6.8070238-10-3 


-9.6117859-10-3 


0.034668583 


-0.082780239 


-0.13230861 


ag 


23.269583 


11.630659 


-0.71127314 


5.3857605 


-0.46984803 


aio 


-0.91120803 


-1.5749174 


4.1212869 


-2.14052 


-0.4187379 


an 


0.2999526 


-0.32863816 


0.12551445 


-0.17836516 


-0.21912773 


an 


-5.3509101 


-2.2770031 


-0.66992889 


-0.64180163 


0.22150197 



Table 9: Fit parameters as defined in Eq. AE.2[l for the functions fqq \Sq5. 





Jqq 


^g^A 




Sqq,3 


oqqA 




1 


1 


-1/3 


1 


1 


a\ 


0.18272533 


0.18271262 


-0.061287765 


0.18253142 


0.18262171 


ai 


0.18635401 


0.18719675 


0.015656292 


0.1756192 


0.18665363 


as 


-0.16337335 


-0.15894226 


-0.07245742 


-0.17432802 


-0.14086689 


a4 


0.26381224 


0.27417865 


-0.089475628 


0.22222852 


0.27777792 


as 


-0.76180191 


-0.78525268 


0.21303639 


-0.6620887 


-0.78622279 


«6 


-0.030248112 


-0.034023724 


0.020746502 


-0.021337505 


-0.033873638 


a-] 


0.012349472 


0.013843671 


-0.022717427 


0.014090792 


0.030078002 


as 


7.7073036-10-3 


7.85052-10-3 


-8.1364483-10-3 


5.9573824-10-3 


8.8545826-10-3 
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Table 10: Fit parameters as defined in Eq. \E.3l for the functions fqg Kg^ 





Jqs 






ai 


0.01532878 


-7.5038708-10-3 


2.1881981-10-3 


a2 


-0.45170812 


-0.29554484 


0.23254702 


as 


-0.1287014 


-0.048548285 


0.022275599 


a4 


0.52370186 


0.30790838 


-0.23459259 


as 


-0.26811216 


-0.15956101 


0.12225428 


ae 


0.022340876 


0.018846548 


-0.017627 


a-] 


-0.19533856 


-0.14177144 


0.11350417 


as 


-0.054689213 


-0.02584021 


0.021535027 


ag 


-0.045358797 


-0.013097499 


1.5210124-10-3 


aio 


-7.7692908-10-4 


1.3753108-10-3 


-1.5369012-10-3 


an 


5.6252092-10-5 


-8.811281-10-5 


9.6359945-10-5 


an 


-8.8107179-10-5 


-5.0153667-10-5 


2.1289443-10-5 


ai-i 


-0.03501338 


-0.013967759 


7.2211643-10-3 



Table 11: Fit parameters as defined in Eq. \E.3l for the functions g 4 3 4 discussed in the text. 













a\ 


1.8239099-10-3 


2.0186389-10-3 


1.6574875-10-3 


1.8413674-10-3 


a2 


0.011452419 


-0.025864066 


-0.25682763 


-0.26685522 


a-i 


-0.015540536 


-0.022761156 


-0.019640631 


-0.025869178 


a4 


-6.5191297-10-3 


0.031656197 


0.25770922 


0.2732 


as 


0.012121019 


-7.5126124-10-3 


-0.13514839 


-0.13361839 


ae 


-4.8820211-10-3 


-3.8773834-10-3 


0.024057235 


0.018832676 


a-] 


-2.9631254-10-3 


-0.021304475 


-0.12327878 


-0.13543713 


a% 


-6.6094881-10-3 


-0.010139352 


-0.018603959 


-0.025245484 


ag 


-3.8199452-10-3 


-4.3853099-10-3 


-4.0261913-10-3 


-6.6003001 - 10 3 


a\o 


1.8172768-10-3 


1.8737765-10-3 


1.065926- 10 3 


9.7391848-10-4 


an 


-1.0969205-10 4 


-1.1280442-10-4 


-5.9556692-10-5 


-5.5377984-10-5 


an 


-2.2922513-10-5 


-2.6895291 - 10-5 


-2.100388-10-5 


-2.4470893 - 10-5 


an 


-2.8989325-10-3 


-5.9168202-10-3 


-4.9177402-10-3 


-6.3797473-10-3 
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kernels, see Eq. (IIV.34I) . In most cases the convolution integrals can be done analytically. For the 
remaining functions we use as ansatz for the process gg ttX: 

4i^3,4 = ai+a2P + a3p2 + a4|33 + a5|3Va6|3ln(|3)+a7|32ln(|3)+a8Pln(P)2 + a9p2ln(|3)2 
+ flioP^ln(p)+aiip3ln(p)2 + ai2pln(p)+ai3p2ln(p)+ai4pln(p)2 

+ai5p^ln(p)2 + ai6ppln(p). (E.4) 

The fitted values for the two functions f . are given in Table fT2l 
For the scale-dependent part in the reaction qg ttX we use 

1^3,4 = 4^,3,4 = P(P'(ailn(p) +a2) +P^a3ln(p) +^4) + p'(a5ln(p) +a6ln(p)') 

+ p(a7ln(p) +a8ln(p)^)) + P% + p^aioln(ri) + p^an ln(ri)2 

+ p3(ai2ln(p) +ai3ln(p)2) +ai4pln(ri) +ai5p3ln(ri) +ai6pln(ri)^ 

+ ai7p^ln(ri)2 + ai8ppln(ri)+ai9p3pln(ri)+a2oPpln(ri)2. (E.5) 

The fitted values for the two functions g[^]^ a are shown in Table [T^ 
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Table 12: Fit parameters as defined in Eq. \E.4l for the two functions g\ 







g^'\ 


a\ 


3 600512- 10""^ 


9 8339551 ■ 10""^ 


a^ 


-0.89103071 


-2.1615978 


as 


4.1432564 


-2.4383759 


a4 


-5.652246 


-1.7863628 


as 


2.396262 


6.3822522 




-0.26139117 


-0.7174869 


av 


0.55643278 


-4.634966 


a% 


-0.029004255 


-0.076074349 


ag 


0.44754783 


0.26832936 


aio 


-0.39236517 


-8.4049829 


an 


-0.52419439 


1.3153541 


«12 


1.9909471 


7.6847524 


«13 


-0.8099637 


-2.9175181 


ai4 


0.013395378 


0.023018999 


«15 


0.031263952 


0.035390606 


ai6 


-1.9746581 


-7.5775663 
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Table 13: Fit parameters as defined in Eq. (IE. 51) for the two functions g\ 





8^'\ 




a\ 


016955545 


-2 7829949- 10-3 




057599714 


03597281 




012017294 


1 5756308 - 10-3 

-1- • «y / «y V_>' ». V V./ K_} J. V./ 


6(4 


-0 061168008 


-0 043921612 


a^ 


025443918 


023078656 


ae. 


-2 4972922- 10"-^ 


-2 5244529- 10-3 




019870969 


019474128 

v_/ « V./ A ^ r / r A \j 


an 


5.0330984- 10^3 


4.949093 - 10-3 


ag 


2.216736-10-3 


7.0180624-10-3 


aio 


7.7028477 - lO^^ 


-1.6586156-10-3 


an 


-1.8013666- 10 4 


1.1083146-10-4 


an 


3.666881-10^ 


-3.7800956-10-^ 


au 


1.7202878-10^3 


1.3536538-10-3 


au 


-4.2292913-10-^ 


2.2933896-10-6 


ai5 


-7.5513607-10-3 


1.7188112-10-3 


ai6 


8.349152-10-^ 


8.6710055-10-^ 


a\i 


1.6975408-10-4 


-1.1585369-10-4 


ai8 


4.1704783-10-*^ 


-2.2181371 - 10-*^ 


a\9 


1.645701-10-3 


-1.7941491 - 10-3 


aio 


-8.3995838-10-^ 


-8.6053362-10-^ 
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